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The question of isospin-symmetry breaking in nuclei of the sd shell is addressed. We propose 
a new global parameterization of the isospin-nonconserving (INC) shell-model Hamiltonian which 
accurately describes experimentally known isobaric mass splittings. The isospin-symmetry violating 
part of the Hamiltonian consists of the Coulomb interaction and effective charge-dependent forces 
of nuclear origin. Particular attention has been paid to the effect of the short-range correlations. 
The behavior of b and c coefficients of the isobaric-mass-multiplet equation (IMME) is explored in 
detail. In particular, a high-precision numerical description of the staggering effect is proposed and 
contribution of the charge-dependent forces to the nuclear pairing is discussed. The Hamiltonian 
is applied to the study of the IMME beyond a quadratic form in the A = ?>2 quintet, as well as 
to calculation of nuclear structure corrections to superallowed O""" — ^ O""" Fermi jS decay, and to 
amplitudes of Fermi transitions to non-analogue states in sd-shell nuclei. 



PACS numbers; 21.60.Cs, 21.10.Pc, 21.10.Jx 

I. INTRODUCTION 

The isospin symmetry is one of the pivotal con- 
cepts in nuclear structure which simplifies largerly many- 
body calculations, for example, within the nuclear shell 
model, and represents a useful guideline in nuclear the- 
ory. The concept is based on the charge- independence 
of the nucleon-nucleon (NN) interaction (invariance un- 
der any rotation in isospin space), which reflects the 
fact that strong proton-proton (w^), neutron- neutron 
{v^n) proton-neutron (w^jj^^x) interactions are to 
a large extent identical. Within the isospin symmetry, 
a many-body nuclear Hamiltonian (without electromag- 
netic interactions) commutes with the isospin operator, 
\H, T] = 0, and its eigenstates can be characterized by an 
isospin quantum number T forming multiplets of (27" -1-1) 
states in a few neighboring nuclei with Tz = — T, . . . , T 
{isobaric analogue states, IAS). 

The charge-independence implies also a charge- 
symmetry of the nucleon-nucleon interaction (invariance 
under rotation by 180° in isospace around the Ty axis), 
. which means the equality of and only, and the 

■ isospin conservation oi N = Z nuclei. This symmetry 
, manifests itself in close similarity of the spectra of mirror 

■ nuclei (up to an overall shift). 

Nevertheless, the isospin symmetry is only an approx- 
imate symmetry in nuclear physics mainly due to the 
Coulomb interaction acting between protons, but also 
due to the presence of charge-dependent forces of nu- 
clear origin. The latter are understood at present to 
have their origin in the difference between the u and d 
quark masses and electromagnetic interactions between 
them [1]. Indeed, the NN scattering data shows unam- 
biguous evidence on the breaking of the two symmetries 
of the NN interaction mentioned above. First, there is a 
small difference between and v^^ (e.g., different scat- 
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tering lengths in the ^5*0 channel: a^ — a^„ = 1.65 ±0.60 
fm [1]) which means the charge-symmetry breaking (or 
charge-asymmetry) of the NN interaction. Second, there 
is an even more substantial difference between and 
on one side and on the other side (e.g., different 
singlet scattering lengths: (a^ -I- a^„) /2 — = 5.6 ± 0.6 
fm) which corroborates the charge-independence breaking 
of the NN interaction. There exist also charge-dependent 
forces which mix the isospin of an NN system, however, 
we do not discuss them here. Detailed consideration 
and theoretical studies of these effects can be found in 
Refs. [1-4]. 

A many-body Hamiltonian containing charge- 
dependent forces does not commute with the isospin 
operator, therefore the isospin is not conserved anymore. 
The eigenstates of such a Hamiltonian represent a mix- 
ture of different isospin eigenstates. This is the case of 
explicit isospin-symmetry breaking and of isospin- mixing 
in nuclear states. 

The degree of isospin nonconservation due to the 
Coulomb interaction and charge-dependent nuclear forces 
is small compared to nuclear effects, however, precise 
description of the isospin-symmetry breaking in nuclear 
states is crucial, when a nucleus is considered as a lab- 
oratory to test the fundamental symmetries underlying 
the Standard Model of the electroweak interaction. One 
of the important applications is the calculation of the 
corrections to nuclear beta decay, which arise due to the 
isospin-mixing in nuclear states and thus should be eval- 
uated within a nuclear many-body model. 

In particular, high-precision theoretical values of nu- 
clear structure corrections to superallowed 0+ — ?> 0+ /3- 
decay rates are of major interest. Combined with various 
radiative corrections, they serve to extract from /i- values 
of these purely Fermi transitions an absolute J^t-vedue of 
the nuclear beta decay. The constancy of ^t for various 
emitters confirms the conserved vector current hypoth- 
esis (CVC) and allows one to deduce the nuclear weak- 
interaction coupling constant, Gp- The ratio of the lat- 
ter with the weak-interaction coupling constant extracted 
from the muon decay gives the absolute value of Vud, 
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the uppcr-lcft matrix element of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix. The upper row of the CKM 
matrix is the one which provides a stringent test for the 
unitarity, while Vud being the major contributor (around 
~94%). The breakdown of the unitarity signifies a pos- 
sibility of new physics beyond the Standard Model, see 
Ref. [5] for a recent review. 

Nowadays, /t-values for thirteen 0+ 0+ /3+ tran- 
sitions among T = 1 analogue states are known with a 
precision better than 0.1%. The largest uncertainty on 
the extracted value (which is of about 0.4%) is due 
to an ambiguous calculation of the nuclear structure cor- 
rection [6]. Therefore, accurate theoretical description of 
isospin mixing in nuclear states is of primary importance. 

Similarly, theoretical calculations of nuclear-structure 
corrections to Fermi /3-decay are necessary to extract the 
absolute .^t value and Vud matrix element from mixed 
Fermi/Gamow- Teller transitions in mirror T = 1/2 nu- 
clei [7]. Nuclear structure corrections to Gamow- Teller 
/3-decay matrix elements arc required in studies of asym- 
metry of Gamow- Teller /3-decay rates of mirror transi- 
tions with the aim to constrain the value of the induced 
tensor term in the axial- vector weak current [8, 9]. 

Apart from the nuclear structure corrections for stud- 
ies of fundamental interactions, precise modclization of 
the Coulomb and charge-dependent nuclear forces is re- 
quired to describe observed mirror energy differences [10] 
and splittings of the isobaric multiplets, amplitudes of ex- 
perimentally measured isospin-forbidden processes, such 
as /3-delayed nucleon emission [11], Fermi /3 decay to non- 
analogue states [12], El transitions in self-conjugate nu- 
clei [13] or isoscalar El component extracted from El 
transitions between analogue states [14] and so on. The 
charge-dependent effective interaction is indispensable for 
understanding the structure of proton-rich nuclei with 
important consequences for astrophysical applications. 

At the same time, accurate theoretical description 
of the isospin-symmetry violation within a microscopic 
model is a great challenge. Various approaches have been 
developed to deal with the problem. 

The first shell-model estimations of isospin mixing are 
dated to the 1960's (e.g. Refs. [15-17]), including their 
applications to the nuclear beta decay studies (see Ref. [8, 
18, 19] and references therein). Among the most recent 
work within the modern shell model, let us refer first to 
the study of Ormand and Brown [20, 21], who constructed 
realistic INC effective Hamiltonians constrained by the 
experimental data (mass splittings of isobaric multiplets). 
Another approach based on the analysis of mirror energy 
differences in p/-shell nuclei was proposed by Zuker and 
collaborators Ref. [22] and gave a profound picture of the 
Coulomb effects. 

It should be remembered that within the shell model, 
one cannot deduce completely a degree of isospin mixing 
in the wave function. The reason is that the Shrodinger 
equation is solved in the harmonic-oscillator basis within 
one or two oscillator shells (valence space) for valence nu- 
cleons only. An INC Hamiltonian allows to introduce the 
isospin-symmetry breaking in the mixing of the many- 
body harmonic oscillator configurations which represent 
Hamiltonian eigenstates. This is sufficient to get the 
energy shifts of isobaric multiplets due to the charge- 
dependent interaction. However, to get matrix elements 



of isospin-forbidden transitions, one has to account for 
the isospin-symmetry breaking beyond the model space. 
With this aim, one has to substitute the harmonic oscilla- 
tor radial wave functions by realistic ones, since the cor- 
rect asymptotics is essential. In this way, the shell model 
allows to predict the rates of isospin-forbidden processes 
which can be compared to experimental data. 

Recent applications of the shell model to superallowed 
/3 decay can be found in Refs. [5, 6, 23, 24] and refer- 
ences therein, while corrections to Gamow- Teller /3 decay 
in mirror systems have been evaluated in Refs. [8, 9]). 
Numerous applications to the isospin-forbidden proton 
emission and to the structure of proton-rich nuclei can 
be found in the literature (e.g. Refs. [25-30]). 

The problem of the isospin-symmetry breaking was in- 
tensively undertaken in the framework of self-consistent 
mean-field theories within the Hartree-Fock -I- Tamm- 
Dankoff or random-phase approximation (RPA) in the 
1990's [31-34]. Recently, more advanced studies have 
been performed within the relativistic RPA approach [35], 
as well as within the angular-momentum-projectcd and 
isospin-projcctcd Hartrec-Fock model [36, 37]. 

Some other many-body techniques have recently been 
applied to deal with isospin non-conservation. In partic- 
ular, evaluation of the isospin mixing in nuclei around 
N Z AO has been performed by variation-after- 
projection techniques on the Hartree-Fock-Bogoliubov 
basis with a realistic two-body force in Ref. [38] . Isospin- 
symmetry violation in light nuclei, applied to the case of 
superallowed decay of the ^°C has been calculated within 
the ab-initio no-core shell model [39], while effects of 
the coupling to the countinuum on the isospin mixing in 
weakly-bound light systems were studied in the Gamow 
shell-model approach [40]. Relation between the isospin 
impurities and the isovector giant monopole resonance 
was explored by Auerbach [41], with a subsequent appli- 
cation to the calculation of nuclear structure corrections 
to superallowed (3 decay [42]. 

Up to now, the approaches mentioned above do not 
agree on the magnitude of isospin impurities in nuclear 
states and predict largely different values for the cor- 
rections to nuclear /3 decay. Given the importance of 
the problematics we have revised the existing INC shell- 
model Hamiltonians. First, since the latest work of 
Ref. [21] there have been accumulated more experimen- 
tal data and data of higher precision on the properties of 
isobaric multiplets (mass excess data and level schemes), 
on isospin-forbidden particle emission, on nuclear radii 
and so on. Development of the computer power and 
shell-model techniques allows us to access larger model 
spaces [43] . In addition, more precise new nuclear Hamil- 
tonians have been designed (e.g. Refs. [44-46]), as well 
as new approaches to accounting for short-range cor- 
relations have been advocated [47, 48]. The purpose 
of this article is to present an updated set of globally- 
parameterized INC Hamiltonians for sd-shell nuclei, and 
to show their quantitative implication to calculations of 
isospin-forbidden processes in nuclei. 

In Section II, we describe the formalism used for a fit 
of the INC interaction. Section III contains the results 
obtained in the sd shell. In section IV we discuss the 
general behavior and numerical agreement of theoretical 
and experimental b and c coefficients, as well as we re- 
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veal and study a so-called staggering phenomenon. In 
Section V wc explore the extension of the IMME beyond 
the quadratic form in the lowest A = 32 quintet. In 
Section VI, we present a new set of nuclear structure cor- 
rections for superallowed 0+ — >■ O"*" Fermi f3 decay, as well 
as a few cases of Fermi transitions to non-analogue states 
(configuration- mixing part). The paper is summarized in 
the last section. 



II. FRAMEWORK 
A. Shell-model Formalism and Fitting Procedure 

Within the nuclear shell model, the eigenproblem is 
solved by diagonalization of the one- plus two-body nu- 
clear Hamiltonian in the basis of many-body Slater deter- 
minants constructed from the single-particle harmonic- 
oscillator wave functions. Since the basis dimension 
grows rapidly with the number of nucleons, the eigen- 
problem is stated only for valence nucleons in a model 
space containing a few (valence) orbitals above a closed- 
shell core. The Hamiltonian matrix thus consists of 
single- particle energies, e^, typically taken from exper- 
iment, and the two-body matrix elements (TBME's) 
Vijki-j- Here indices i, j, fc, I denote full sets of quantum 
numbers (nlj) necessary to characterise a given single- 
particle orbital, while J denotes the total angular mo- 
mentum of a coupled two-body state. 

Wc suppose that proton-proton, neutron-neutron and 
proton- neutron matrix elements may all be different. 
Similarly, proton and neutron single-particle energies are 
not the same. The goal is to find an interaction which 
describes well both nuclear structure and the splitting 
of isobaric multiplets of states. In principle, an effective 
shell-model interaction may be derived microscopically 
from the bare NN force by applying a renormalization 
technique [49, 50]. However, such interactions, obtained 
from a two-body potential only, should still be adjusted, 
in particular, to get correct monopole properties [43, 51]. 
This is done by a least-squares fit of the monopole part of 
the Hamiltonian or of the whole set of TBME's to exper- 
imental data. Since the number of the matrix elements is 
huge, it is not feasible for the moment to get a realistic 
charge-dependent effective interaction in this way. 

An alternative approach to the problem is first to get 
a reliable effective shell-model interaction in the isospin- 
symmetric formalism adjusted to describe experimental 
ground and excited-states energies. Then to add a small 
charge-dependent part within perturbation theory and 
to constrain its parameters to experimental data. Diago- 
nalization of the total INC Hamiltonian in the harmonic 
oscillator basis will lead to isospin mixing. 

In the sd shell- model space (consisting of the 0^5/2, 
lsi/2, and 0(1^/2 orbitals) the most precise isospin- 
conserving Hamiltonians, denoted below as H, are the 
USD interaction [52], as well as its two more recent ver- 
sions USDA and USDB [44]. First, we obtain its eigen- 
values and eigenvectors: 

H\a,T,T,) = (Ho + Vo)\a,T,T,) = E{a,T)\a,T,T,) . 

Here, a = {A, J"^ , N^xc, ■ ■ ■) denotes all other quan- 
tum numbers (except for T and T^), which are re- 



quired to label a quantum state of an isobaric multiplet. 
E{a,T) is independent from Tr. iJo is the independent- 
particle harmonic oscillator Hamiltonian which involves 

the (isoscalar) single-particle energies el"^ = (ef -f £")/2, 
while Vq stands for a two-body residual interaction in the 
sd shell. 

Then we construct a realistic isospin-symmetry violat- 
ing term to get a total INC Hamiltonian. In general, we 
consider a charge-dependent interaction, which includes 
the Coulomb interaction acting between (valence) pro- 
tons, and also charge-dependent forces of nuclear origin. 
The Coulomb interaction reads 



Vcoui{r) 



(1) 



while the charge-dependent nuclear forces are represented 
in this work either by a set of scaled T = 1 matrix ele- 
ments of the isospin-conserving interaction Vb (denoted 
as Vq^^^) or by a linear combination of Yukawa- type po- 
tentials: 



Vf,{r) 



exp (/i^r) 
exp (/ipr) 



(2) 



where /^tt — 0.7 fm^^ and jip = 3.9 fm~^, correspond- 
ing to the exchange of pion or p meson, respectively, 
and r being the relative distance between two interact- 
ing nucleons. The Coulomb interaction contributes only 
to the proton-proton matrix elements, while the charge- 
dependent nuclear forces may contribute to all nucleon- 
nucleon channels. Thus, we can express the charge- 
dependent part of the two-body interaction as 



= ^couiycoul{'r) 



q—pp.mi.p7i{T—l) 



{\lV:^{r)+\lV;}{r)+\lVS) (3) 



where Vq is the same as Vq^^, while XcouU ^1, A^, Aq 
are strength parameters characterizing the contribution 
of charge-dependent forces. These parameters can be es- 
tablished by a fit to experimental data. 

The two-body charge-dependent interaction V in 
Eq.(3) can alternatively be decomposed in terms of ranks 
0, 1, and 2 tensors in the isospin space as 

The corresponding two-body matrix elements can be re- 
lated to those in proton-neutron formalism, i.e. 



(0) 



ijkl 



(1) 



ijkl 



(2) 



ijkl 



yPP 



yPP 
yPP 



,/ + ^ijkl 
^ijkl-J ' 
^ijkl-J 



y 



pn{T=l) 



ijkl: J 



(4) 



2V: 



pn(T=l) 
ijkl; J 



In addition, the charge-dependent part of the Hamilto- 
nian may contain a one-body term, H^j^ of a pure isovec- 
tor character, which involves the isovector single-particle 



energies (ISPE's), e. 



(1) 



This term accounts for 
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the Coulomb effects in the core nucleus. Thus the most- 
general charge-dependent part of the effective Hamilto- 
nian reads 



He 



D 



H, 



16 
CD 



V . 



The charge-dependent part of the effective interaction 
is well known to be small and to be mainly of two-body 
type. The shift of isobaric multiplets due to the presence 
of charge-dependent Hamiltonian, Hcd, in lowest order 
of perturbation theory is given by its expectation value in 
the states having good isospin: (a,T,Tz \ Hcd \ct,T,Tz). 
Application of the Wigner-Eckart theorem leads to the 
following expression: 

{a,T,T^\ Hcd \a,T,T^) 
+E^^\a,T)T, 

+E^^^ {a, T) [3T2 - T(r + 1)] , (5) 

where the isoscalar part T^^"-* contributes only to the over- 
all shifts of the multiplet, the isovector part F'^-'^' and 



ISPE's (e^i)) resuhs in E^'^\a,T), while the isotensor 
part is the only contributor to E'^'^^a.T). The lat- 
ter two terms lead to the splitting of the isobaric multiplet 
and to the isospin mixing in the states. 

Based on this asumption, Wigner showed [53] that a 
quadratic isobaric-mass-multiplet equation (IMME), 

M{a, T, T,) = a(a, T) + 6(a, T)T, + c{a, T)T^ , (6) 

is sufficient to approximate the splitting of isobaric mass 
multiplets for a given a, and T. The a, b and c are 
coefficients. 

Since only the isovector and isotensor part of Hcd 
could lead to isospin-symmetry violation (to splitting of 
the isobaric multiplets and to isospin mixing) , we will be 
interested in these two terms only. Futhermore, in the fit 
of the nuclear TBME's in the isospin-symmetric formal- 
ism, part of the isoscalar Coulomb term has been taken 
into account by an empirical correction to the experi- 
mental binding energies (see [44] and references therein) . 
Therefore, we add to the isospin-conscrving Hamilto- 
nian a charge-dependent Hamiltonian, containing isovec- 
tor (iv) and isotensor (it) terms only, namely. 



9=1,2 



q=l,2 V 
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CD 



(7) 



where q now denotes the isotensor rank of the operators 
and labels the corresponding strength parameter, while 
the indice v is used to list all separate terms. The second 
line of Eq.(7) includes the one-body term, with being 
the corresponding strength parameters. 

The isovector £^(^^(a, T) and isotensor E'-^^{a,T) con- 
tributions to the expectation value of H^^^* (or Hcd), 
can be either extracted from the energy shift due to the 
isovector F^^) (or H^^jj) and isotensor F'^' parts of the 
charge-dependent Hamiltonian, respectively, or from cal- 
culations of the energy shifts of all multiplet members. 
Following the latter method, we represent the TBME's 
of Vi, in terms of the proton-proton matrix elements only 
and then we calculate its expectation value in each state 
of the multiplet E^{a,T,T^) = (a, T, T^l K |a, T, T,). 
Then, the isovector and isotensor contributions to a given 
multiplet of states are respectively expressed as 



i?W(a,r) 
El'\a,T) 



T{T+1){2T 



— J2 {-T.)EAc^,T,T,) , 



T^ = -T 



T{T + 1)(2T - 1){2T + 1)(2T -t- 3) 

T 

X J2 [iT^^T{T+l)]E,{a,T,T,). (8) 

The same method holds also for the ISPE's. Summing 
over all contributions to Hcd, we get theoretical IMME 



b and c coefficients as 

b"'ia,T) = j:xi'^Ei'\a,T), 

c*''{a,T)^3Y.xi^^Ei^\a,T) 



(9) 



ISPE's are only included into the expression for b val- 
ues. 

To find the best strengths we have performed a 
least-squares fit of theoretical 6*'' and c*'' coefficients to 
experimental IMME b and c coefficients : b^^^ i ai {i = 
1, . . . , Nb) and c^j^^ ± aj (j = 1, . . . , Nc). Implying that 
they have a Gaussian distribution, we have minimized the 
deviation (e.g., for b coefficients): 



E 



exp 



hth\2 



with respect to the parameters xl}\ i.e. 



dX^ 



d 



(1) 



dx\}^ It 



exp 



Uth\2 



(10) 



(11) 



which has lead us to a system of linear equations for Xi, 

N, 



(1) 



E 



E 



\(i)p(i)p(i) 



= 0. 



In matrix form this system looks as 
AW = Q or ^ KW, 



(12) 



(13) 
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FIG. 1. (Color online) Harmonic oscillator energy spacing, hw: (i) Blomqvist-Molinari formula [54] (solid black line), (ii) 
Blomqvist-Molinari formula with an additional scaling factor given by Eq.(3.7) of Ref. [21] (double-dot-dashed blue line); (iii) 
Blomqvist-Molinari formula refitted by Kirson [55] (dashed red line); (iv) huj deduced from the experimental nuclear charge 
radii [56] following the procedure indicated in Ref. [55] with fi^/m — 41.458 MeV.fm^ (black dots). Only A — 2, . . . , 60 is 
indicated. 



with 

1=1 

Nb 771(1) ^exp 

Q.-T.^^- (14) 

i=l ' 

Since theoretical b and c coefficients are linear functions 
of the unknown parameters A^'^ of Eq.(9), the fitting pro- 
cedure is reduced to solving linear equations. Solution of 
these equations with respect to A results in a set of the 
most optimal strength parameters aI.^^: 

A = QW-^ (15) 

To get uncertainties of the strength parameters found, 
we evaluate the root-mean-square (rms) deviation from 
the error matrix as 

aaW = \/((aW -aF)2) - ^{w-^U ■ (16) 

A similar procedure holds for the adjustment of c coeffi- 
cients. 



After adjusting the interaction, we solve the eigen- 
problem for a thus constructed INC Hamiltonian in the 
proton-neutron formalism: [Hi^c ^T] ^ 0: 

HiNc\oip,a.n) = (iJ + -ffg'^''*)|ap,a„) = E\ap,an) ■ 

As a result, the Hamiltonian eigenstates do not possess 
good isospin quantum number anymore and thus arc mix- 
tures of different T values. 



B. TBME's of the Coulomb and Yukawa- Type 
Potentials 

1. Harmonic Oscillator Parameter 

The TBME's of the Coulomb and Yukawa-type po- 
tentials Eq.(l) and Eq.(2), used to calculate the energy 
shifts, were evaluated using the harmonic-oscillator wave 
functions for mass A — i% and the subsequent scaling 

In Ref. [21], hw was taken in its most commonly used pa- 
rameterization expressed by the Blomqvist-Molinari for- 
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FIG. 2. (Color online) Vcouiir) without SRC (dotted black line), Vcoui(f') with proposed parameters for Jastrow type SRC 
on the basis of coupled-cluster calculation with CD-Bonn (double-dashed red line) or with Argonne V18 [48] (solid blue line), 
Vcouiij) with Miller-Spencer parameterized Jastrow type SRC [57] (dashed purple line); Vcoui(^+('')) with UCOM SRC (dot- 
dashed green line), Vcoui(i?+'('')) with UCOM SRC (double-dot-dashed brown line). The inset enlarges the left-hand part of 
this figure. 
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-1/3 



25^ 



-2/3 



(18) 



For the sd shell, an additional scaling factor was imposed 
(see Eq. (3.7) in Ref. [21]) to improve the agreement with 
the data at the beginning and at the end of the sd shell. 

However, recent empirical values of hw. derived from 
updated experimental nuclear charge radii in Ref. [56], 
differ significantly from the values predicted by Ormand 
and Brown in Ref. [21], especially in the middle of the sd 
shell, not considered in the latter work. The comparison 
is shown in Fig. 1. Some improvement is reached by a 
recent global parameterization of the Blomqvist-Molinari 
formula for the whole nuclear chart (A = 2, . . . , 248) per- 
formed by Kirson [55] (see Fig. 1). 

We have performed a fit with both parameterizations of 
hiij, however, none of them resulted in a sufficiently low 
rms deviation values in our fit for h and c coefficients. 
The plausible reason is that the existing parameteriza- 
tions for hw values in the sd shell are not close to the 
values extracted from experimental nuclear charge radii. 
To overcome this difficulty, in the present work we have 
scaled the TBME's as given by Eq.(17), directly using 
experimentally based values for hijj values in the sd shell, 
mentioned above and shown in Fig. 1. The ISPE's were 



also evaluated for A 
Eq.(17). 



39 and then scaled as given by 



2. Short-Range Correlations 

Since the TBME's of Coulomb or meson-exchange po- 
tentials are calculated by using harmonic-oscillator wave 
functions, it is important to account for the presence of 
short-range correlations (SRC). We have carefully studied 
this issue by two different methods. First, the Jastrow- 
type correlation function, which modifies the relative part 
of the harmonic-oscillator basis, (r), to 

<P'ni{r) = [l + f{r)] Mr), 

with /(r) being parametrized as 

/(r) = -7e-"'-'(l-/3r2) . (19) 

Then the radial part of the TBME's of the Coulomb and 
Yukawa type potentials between the modified harmonic- 
oscillator wave functions (f>ni{r) and ^^^/((f) becomes 

oo 



(20) 
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We used three different sets of parameters a, /? and 7 in 
Eq.(19): those given by Miller and Spencer [57] and two 
alternative sets recently proposed on the basis of coupled- 
cluster studies with Argonne (AV18) and CD-Bonn po- 
tentials [48] (see Tabic I). For brevity, we will refer to the 
two latter sets as CD-Bonn and AVIS. 



TABLE I. Parameters for the short-range correlation function. 





a 


P 


7 


Miller-Spencer 


1.10 


0.68 


1.00 


CD-Bonn 


1.52 


1.88 


0.46 


Argonne-V18 


1.59 


1.45 


0.92 



It is seen that the Miller-Spencer approach to SRC 
quenches the Coulomb matrix element (as well as that 
of the Yukawa p-meson exchange potential) and thus 
reduces Coulomb expectation values more compared to 
other SRC schemes. Interestingly, the CD-Bonn parame- 
terization and in some cases the AV-18 parameterization 
show even a small increase of the Coulomb expectation 
value. Bearing this in mind, in the next section we, how- 
ever, perform a fit of the INC parameters for all cited 
approaches to the SRCs. 

III. RESULTS AND DISCUSSION 



Besides, we have also used another renormalization 
scheme following the unitary correlation operator method 
(UCOM) [58]. Since we need to correct only central op- 
erators, the UCOM reduces to the application of central 
correlators only, i.e. the radial matrix elements are of the 
form 



(i>n'i{r)v {R+{r)) 4>ni{r)dr , 



(21) 



where two different R+{r) functions have been used in 
S = Q,T = 1 and 5 = 1, T = 1 channels, namely. 



R\{r) = r + a\-] exp 



-expij 



(22) 



with a = 1.3793 fm, jS ^ 0.8853 fm, ry = 0.3724 in the 
S" = 0, T = 1 channel, and 



i?"(r) 



a [ 1 — cxp ( 



exp 



1.3888 fm, 7 



oxplj^ 

(23) 

0.1786 in the 



with a = 0.5665 fm, /3 
5 = l,r= 1 channel [58]. 

The modifications of Vcoui brought about by different 
approaches to the SRC issue are shown in Fig. 2. Similar 
trends hold for Vp and V^. 

Although the UCOM renormalization scheme differs 
from the Jastrow-type correlation functions, we can easily 
notice that either of the i?+ (r) functions does not strongly 
affect the original potentials. Somewhat stronger modi- 
fications are brought about by the CD-Bonn based pa- 
rameterization. The Miller-Spencer parameterization of 
the correlation function induces the highest suppression 
of the potentials at short distances and leads to a van- 
ishing value at r = 0. Similar conclusions are reported 
in Ref. [59] in the context of double-beta decay studies. 
Strong modifications are clearly seen for AV18 as well. 

To illustrate the effect from different approaches to the 
SRCs on the results to be discussed in the following, we 
present in Table II the ratios of the Coulomb expecta- 
tion values in the ground and several low-lying excited 
states of a few selected nuclei from the bottom, from the 
top, and from the middle of the sd shell-model space, i.e. 
^^Ne (2 valence protons), ^^K (2 proton holes), and ^°S 
and ^^Mg, respectively. The second column of Table II 
contains absolute expectation values of the bare Coulomb 
interaction, while the other columns show the ratios to 
it from Coulomb interaction expectation values which in- 
clude SRC. 



A. Fitting Procedure 

We have followed the fitting strategy proposed in 
Ref. [21]. First, we construct theoretical b and c coef- 
ficients Eq.(9) as described in the previous section (using 
experimentally based hu and accounting for the SRC by 
one of the above mentioned methods). Then, we sep- 
arately fit them to experimental b and c coefficients to 
get the most optimal values of aI.^"* and respectively. 
We assume here that the isovector and isotcnsor Coulomb 
strengths arc equal. To this end, the isovector and isotcn- 
sor Coulomb strengths obtained in both fits are averaged 



(A 



coul 



(A 



(1) 



(2") 

Aco«;)/^) ^'^'^ ^'''^ kept constant. Then 



the rest of the strength parameters are refitted with this 
fixed Coulomb strength. 

In order to verify our method, we performed a direct 
comparison with the results of Ref. [21]. We have fol- 
lowed their setting exactly by adopting the experimental 
values from Table 5 [60] of Ref. [21], the parameterization 
of the hio and the scaling factors (see Eqs.(3.5)-(3.7) of 
Ref. [21]) for TBME's of V and ISPE's, as well as the 
Miller and Spencer Jastrow-type function [57] to account 
for the SRC effects [61]. We have also imposed certain 
truncations on calculations for ^ = 22 and A — 34, as 
was done in that work [21]. 

In this way, we have successfully reproduced the 
strength parameters given in Table 2 of Ref. [21]. 

For curiosity, besides the USD interaction, we have also 
tested USDA and USDB [44] , keeping the number of data 
points as selected by Ormand and Brown, but using up- 
dated experimental values from Ref. [62]. No truncations 
were used in the calculations for ^ 22 and A ~ 34. The 
corresponding strength parameters are given in Table III. 
The uncertainties on the strength parameters have been 
deduced from Eq.(16). They are significantly smaller 
than the values published in Ref. [21] due to the fact that 
the authors used some folding with the rms deviation [63] . 
It is remarkable that there is no much difference between 
various nuclear interactions for the small data set and 
all parameter strengths are in agreement with the range 
of values found by Ormand and Brown (uncertainties in- 
cluded). 



B. Experimental data base of b and c coefficients 

In the present study we use for the fit an extended and 
updated experimental data base where all latest relevant 
experimental mass measurements and excited states have 
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TABLE II. Ratios of Coulomb expectation values of ^*Ne, '^*K, '^"S, ^^Mg, produced from various SRC approaches to evaluation 
without (w/o) SRC. 







(''i/'l 1/ 7 / e-Do 1?/;^ 

^coulw/o bHC |y^/ 

(MeV) 


Miller-Spencer 


{ipl Vcoul vjith SRC 

CD-Bonn 


IV'} / (V'l Vcoul^/o SRC 1^) 

Argonne V18 


UCOM 


Mass 18, 


i»Ne 
















0"'' e- s 




0.531 


0.900 


1 


.008 




0.978 


0.958 


2+ 




0.449 


0.961 


1 


.010 




0.997 


0.981 


4+ 




0.389 


0.984 


1 


.007 




1.001 


0.991 


0+ 




0.412 


0.952 


1 


.005 




0.990 


0.979 


2+ 




0.380 


0.993 


1 


.006 




1.003 


0.994 


0+ 




0.425 


0.980 


1 


.011 




1.004 


0.987 


Mass 38, 


















0+ E s 




16.402 


0.986 


1 


.007 




1.003 


0.991 


2+ 




16.316 


0.986 


1, 


.007 




1.003 


0.992 


Mass 30, 


30g 
















U g.b. 




10.721 


0.984 


1, 


.007 




1.002 


0.990 


2+ 




10.696 


0.985 


1, 


.007 




1.002 


0.991 


2+ 




10.704 


0.985 


1, 


.007 




1.002 


0.991 


1+ 




10.632 


0.987 


1, 


.007 




1.003 


0.992 


Mass 26, 


2«Mg 
















0+ g.s. 




2.518 


0.967 


1 


.008 




0.998 


0.984 


2+ 




2.480 


0.974 


1 


.008 




1.000 


0.986 


2+ 




2.491 


0.972 


1, 


.008 




0.999 


0.986 


0+ 




2.491 


0.972 


1 


.008 




0.999 


0.986 



TABLE III. Strength parameters obtained in a fit to the num- 
ber of data points selected as in Ref. [21]. 



USD 



USDA 



USDB 



rms (keV): 
b coefficients 
c coefficients 

-A[,^^ X 100 
-A^^' X 100 
4d5/2 (MeV) 
4«/2 (MeV) 



-031/2 



(MeV) 



23.3 
6.9 

1.0077 (1) 
1.3430 (53) 
4.0473 (122) 
3.4076 (2) 
3.3269 (6) 
3.2739 (4) 



28.7 
8.4 

1.0157 (2) 
1.2284 (64) 
5.1755 (146) 
3.4062 (2) 
3.2966 (6) 
3.2853 (5) 



26.8 
8.8 

1.0168 (2) 
1.4669 (64) 
4.9180 (139) 
3.4009 (2) 
3.2898 (6) 
3.2756 (5) 



in a full sd shell-model space. 

• Range I. It includes all ground states (g.s.) and a 
few low-lying excited states throughout the sd shell 
(note that the middle of sd shell was not considered 
in Ref. [21]). This range consists of 81 6 cocfEcients 
and 51 c coefficients. For excited states, the discrep- 
ancy between the energy calculated by the isospin- 
symmetry invariant Hamiltonians and experimental 
excitation energy is less than ~200 keV. 

• Range II. It represents an extension of Range I, 
which includes more excited states. It contains 26 
more T = 1/2 doublets, an additional triplet and 
an additional quartet of state resulting in 107 b co- 
efHcients and 53 c coefficients. 



been taken into account. Indeed, in Ref. [21], the selected 
experimental data consisted of the bottom (A = 18 — 
22) and the top {A = M - 39) of the sd shell-model 
space and included 42 experimental b coefficients and 26 
experimental c cocfficents. 

To get a realistic INC Hamiltonian, we take into ac- 
count in the present fit all available and well-described by 
the sd-shell model isobaric doublets (T = 1/2), triplets 
(T = 1), quartets {T = 3/2) and quintets (T = 2) for 
nuclei between A = 18 and A = 39. The experimentally 
deduced values of the IMME a, 6, c (and d, e) coefficients 
are taken from Ref. [62], which represent an up-to-date 
version of the previous evaluation performed by Britz et 
al. [64]. In particular, the revised experimental database 
incorporates results of all recent mass measurements from 
the evaluation [65] (or given in specific references) and 
modern experimental level schemes [66]. 

In this work we have used three different ranges of data 



• Range III. The widest range, which tops up Range 
II with 32 more excited states from 25 doublets, 6 
triplets, and an additional quartet, resulting alto- 
gether in 139 b cocfficents and 60 c cocfficents. 

These three ranges of selected experimental data points 
are the same for each fit with either the USD, or USDA, 
or USDB interactions. They are presented and discussed 
in Section IV. 



C. Results of the Fit 

All calculations have been performed in an untrun- 
cated sd shell. The TBME's of the schematic inter- 
actions (Coulomb and meson exchange potentials) have 
been evaluated for ^ = 39 and scaled using experimen- 
tally obtained huj. The fit procedure is stated in Sec- 
tion III A. 
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1. INC Hamiltoman and Coulomb strength 

We have tested five different combinations of the effec- 
tive charge-dependent forces: (i) Vcouh (ii) Vcoui and V^, 
(iii) Vcoui and Vp, (iv) Vcoui and Vq, (v) Vcoui, Vp, and 14 ■ 
The main criterium for the choice of the best Hamilto- 
nian structure was the value of the rms and the value of 
the Coulomb strength which was kept as a free parame- 
ter. It turned out that almost all combinations gave sim- 
ilar rms values (within 2 keV). However, on the basis of 
the Coulomb strength parameter we could make a selec- 
tion. We suppose here that the Coulomb strength should 
be close to unity. Indeed, higher-order Coulomb effects 
which are not taken into account here may be respon- 
sible for some deviations of the Coulomb strength from 
unity. However, we suppose that this may be within 1-2% 
and any stronger renormalization (5% or more) should be 
avoided. 

The Coulomb strengths from various combinations of 
the INC Hamiltonians are summarized in Fig. 3. The 
calculations correspond to the USD interaction and a fit 
to the data from Range I, while all approaches to the SRC 
were taken into account. Other choices of the isospin- 
conscrving interaction and other ranges of data selections 
produce similar trends and results. 

First, using the Coulomb interaction as the only source 
of the isospin-symmetry breaking produces a reasonable 
value of the isovector strength (around 1.00), but the 
isotensor strength largerly deviates from unity (up to 
1.19), with the corresponding rms deviations of around 
36 keV for b coefficients and of around 18 keV for c co- 
efficients. The average Coulomb strength, Xcoui is there- 
fore larger than unity (around 1.10) and results in an 
increased rms deviation for b coefficients. The resulting 
parameter strength are summarized in Table IV. This is 
the manifestation of the so-called Nolen-Schiffer anomaly 
first evidenced in T = 1/2 mirror energy shifts [67] and 
later also found in T = 1 displacement energies (e.g. 
see Ref. [68, 69] and references therein). We find that 
the Coulomb potential alone satisfactorily describes the 
mirror energy differences (low rms deviations for b co- 
efficients), possibly due to the fact that the Coulomb 
effects of the core are taken into account through em- 
pirical ISPE's (established by the fit as well). However, 
the Coulomb force alone does not reproduce experimen- 
tal isotensor shifts (larger values of the Coulomb isoten- 
sor strength). Since the sd shell- model wave functions 
include configuration mixing fully within the Ohuj model 
space, this may be an evidence for the necessity of charge- 
dependent forces of nuclear origin. 

Next, it turns out that the Coulomb interaction com- 
bined with the pion-exchange potential Vt^ also requires 
a strong renormalization of the Coulomb strength. This 
was noticed already by Ormand and Brown in Ref. [21]. 
The Coulomb strength reduces to about 0.8 for the Miller- 
Spencer parameterization of the Jastrow function, while 
this factor is around 0.9 — 0.95 for other SRC approaches. 
For these reasons, we do not use pion exchange to model 
charge-dependent nuclear forces in this work. 

A better description is provided by the exchange of a 
more massive meson, e.g. the p meson. Following the- 
oretical studies [70-72], we use in the present work an 
85% reduction in the mass of p meson. A better agree- 



ment with the exchange of a meson heavier than the pion 
may signify a shorter range of a charge-dependent force 
of nuclear origin. 



colli 



1 

2 
3 
4 

3 
4 

2 
3 
4 
5 



0.80 0.85 0.90 0.95 1.00 1.05 _ 1.10 
Averaged Coulomb Strength Parameters X^^^^j 

FIG. 3. (Color online) Average Coulomb strength parame- 
ter, Acoui, as obtained from the fit with the USD interaction 
to the Range / data selection (see Section III A for details). 
Down (blue) triangles correspond to the fit with Coulomb 
force alone. The Xcoui obtained from Vcoui and 14 are de- 
picted by (purple) dots, and up (black) triangles are \coui ob- 
tained from Vcoui and Vp combination, whereas (red) squares 
represent Xcoui from a fit with the Vcoui and Vb combination 
of the two-body charge-dependent forces. Y-axis tic labels: 1, 
"without SRC"; 2, Miller-Spencer; 3, CD-Bonn; 4, AV-18; 5, 
UCOM. 

We confirm also the conclusion of Ref. [21] that a com- 
bination of the pion and p meson exchange potential to 
model nuclear charge-dependent forces does not allow one 
to improve the value of the rms deviation. This is why we 
present here strength parameters only for two combina- 
tions of the charge-dependent forces from the list above, 
namely, (iii) Vcoui and Vp and (iv) Vcoui and Vq. The re- 
sulting rms deviation of these fits and the corresponding 
Coulomb strengths are indeed rather close, in agreement 
with the conclusion of Ref. [21]. We discuss both cases in 
the next section. 



2. rms Deviation Values and Strength Parameters 

Table V gives an overview of strength parameters for 
two types of the INC Hamiltonian: (iii) Vcoui and Vp 
(columns 3, 5, and 7) and (iv) Vcoui and Vq (columns 2, 
4, and 6). Calculations have been performed with the 
USD, USDA, and USDB nuclear Hamiltonians and for 
each of the three data ranges. All four approaches to SRC 
(Jastrow type function with three different parameteriza- 
tions or UCOM) from Section II B 2 have been tested and 
the intervals of parameter variations are indicated in the 
table. 

As seen from Table V, the rms deviation changes little 
for various types of the SRC (within 1 keV) and for both 
types of the charge-dependent Hamiltonian. 

The rms deviation turns out to depend mainly on the 
number of data points used in a fit. It is remarkable 
that although Range I contains almost twice the num- 
ber of data points of Ref. [21], the rms deviation increases 
only by ~5 keV. Overall, the rms deviation of Range 
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TABLE IV. Fitted 


strength parameters" 


for Coulomb as an only 


source of the isospin-symmetry breaking 


; force. 




w/o SRC 


Miller-Spencer 


USD 
CD-Bonn 


Argonne V18 


UCOM 


rms (keV): b coefficients'' 
rms (keV): c coefficients 

41/2 (MeV) 

s^^,,, (McV) 


60.4 

20.1 

1.074 

3.234 

3.060 

3.261 


79.9 

26.5 

1.110 

3.216 

2.963 

3.265 


DU.O 
20.0 

1.067 
3.234 
3.057 

3.262 


DO. 8 

21.8 
1.078 
3.229 
3.030 

3.263 


67.6 

22.5 

1.091 

3.228 

3.025 

3.263 



^ All strength parameters are presented in 4 significant figures. 

^ Before averaging A^^^j and a'^^^j, tlie rms deviations for b coefficients are ~32.8 keV. 



TABLE V. Various combinations of INC potential and their strength parameters. 



Ran (TP 

J-V dud XLCtXi.^^ 


USD 


USDA 


USDB 




Vcoul + Vo 


Vcoul + Vp 


Vcoul + Vo 




Vcoul + Vo 


Vcoul + Vp 


Range I 
rms (keV): 
b coefficients 


-32.2 


-32.7 


-35.2 


-34.5 


-34.3 


-33.9 


r'Opflrir'ip'n'J'C! 


~9.1 


-10.5 


-10.3 


-10.9 


-9.7 


-10.6 




1.006 - 1.015 


1.002 - 0.9839 


1.007 1.016 


0.9993 - 0.9852 


1.008 - 1.018 


0.9987 - 0.9845 


-A^^' X 100 


0.8162 - 1.637 




0.7845 - 1.715 




0.9211 - 1.815 




-A^^^ X 100 


2.629 - 3.861 




2.915 - 4.227 




2.851 - 4.122 




-A^^^ (MoV) 


— 


4.465 - 100.0 


— 


0.8517 - 82.69 


— 


2.228 - 89.52 


A^^^ (McV) 


— 


33.52 - 209.8 


— 


35.27 - 217.1 


— 


35.40 - 216.9 


^11/2 (MeV) 


3.278 - 3.279 


3.294 - 3.295 


3.269 - 3.276 


3.295 - 3.298 


3.267 - 3.272 


3.293 - 3.295 


e</j,,, (MeV) 


3.277 - 3.299 


3.294 - 3.302 


3.273 - 3.298 


3.301 - 3.312 


3.265 - 3.286 


3.297 - 3.306 


e*^' (MeV) 


3.319 - 3.336 


3.344 - 3.346 


3.327 - 3.346 


3.360 - 3.367 


3.323 - 3.341 


3.358 - 3.362 


Range II 
rms (keV): 
b coefficients 


-44.1 


-44.1 


-46.4 


-47.0 


-45.5 


-46.6 


r'Opfnr'ipTi^c; 


-9.3 


-10.6 


-10.4 


-10.9 


-9.8 


-10.7 


^coul 


1.008 - 1.017 


0.9808 - 1.005 


1.008 - 1.017 


0.9826 - 1.006 


1.009 - 1.019 


0.9814 - 1.005 


-X^o^ X 100 


1.202 - 2.019 




1.149 - 2.079 




1.324 - 2.215 




-A^^^ X 100 


2.611 - 3.843 




2.901 - 4.213 




2.836 - 4.106 




-A^^^ (MeV) 




10.44 - 120.3 




6.601 - 103.3 




8.366 - 110.8 


A^^^ (McV) 




33.87 - 212.2 




35.54 - 219.1 




35.76 219.3 


4%/2 (MeV) 


3.271 - 3.273 


3.289 - 3.291 


3.267 - 3.271 


3.292 - 3.295 


3.257 - 3.261 


3.261 - 3.265 


4%/2 (MeV) 


3.273 - 3.295 


3.283 - 3.290 


3.271 - 3.296 


3.289 - 3.299 


3.266 - 3.288 


3.262 - 3.282 


4ti/2 (MeV) 


3.283 - 3.301 


3.300 - 3.304 


3.290 - 3.310 


3.313 - 3.322 


3.320 - 3.342 


3.286 - 3.305 


Range III 
rms (keV): 
h coefficients 


-65.0 


-65.2 


-67.4 


-67.4 


-65.7 


-66.0 


c coefficients 


-10.2 


-10.5 


-11.3 


-10.8 


-10.6 


-10.5 


Xcoul 


1.015 - 1.025 


0.9743 - 1.004 


1.017 - 1.027 


0.9780 - 1.007 


1.018 - 1.028 


0.9752 - 1.005 


-A^^' X 100 


2.331 - 3.189 




2.453 - 3.423 




2.591 - 3.515 




-A^^' X 100 


2.408 - 3.648 




2.690 - 4.012 




2.634 - 3.915 




-A^^^ (MeV) 




37.85 - 228.6 




3.300 - 210.6 




36.02 - 220.4 


A^^^ (McV) 




33.61 215.0 




35.04 - 221.1 




35.54 - 222.5 


41/2 (MeV) 


3.245 - 3.247 


3.256 - 3.258 


3.237 - 3.241 


3.261 - 3.262 


3.233 - 3.237 


3.260 - 3.261 


4%/2 (MeV) 


3.228 - 3.251 


3.169 - 3.186 


3.226 - 3.251 


3.179 - 3.191 


3.219 - 3.240 


3.174 - 3.187 


4ti/2 (MeV) 


3.152 - 3.168 


3.125 - 3.127 


3.147 - 3.165 


3.125 - 3.127 


3.153-3.170 


3.127 - 3.131 
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// is ~30% higher compared to Range I, while the rms 
deviation value for Range III is about twice as large 
as that of Range I. It should also be remembered that 
low-lying states calculated with the isospin-conserving 
USD/USDA/USDB interactions arc in general in better 
agreement with experiment than high-lying states. 

We notice that the USD interaction always produces 
slightly lower rms deviations than USDB and USDA. This 
happens even in the fits to Range III data, although the 
USD was adjusted to a smaller set of excited levels as 
compared to the later versions USDA and USDB. 

Variations in the values of the parameters indicated in 
each entry of the table are due to the different SRC ap- 
proaches. In general, more quenched expectation value of 
an operator results in a higher value of the corresponding 
parameter strength. The most crucial role is played by 
the Coulomb potential, since it is the major contribution 
to isobaric mass splittings. Deviations can be slightly 
greater or less than unity for different combinations of 
charge-dependent forces. 

To reduce the discrepancy, the strengths of the charge- 
dependent forces of nuclear origin, x'^^^^^.^i, are adjusted 
in the fit in a way to match experimental isobaric mass 
splittings. We keep the isovcctor and isotcnsor strengths 
of the nuclear charge-dependent forces as two indepen- 
dent parameters. 

fl- Vcoui o-nd Vo combination. These combinations 
almost always produce the lowest rms deviations for b and 
c coefficients. Fitted to the smallest range of data, the 
isovector and isotensor strengths of the nuclear isospin- 
violating contribution represent about 0.7-1.7% and 2.9 
- 4.2%, respectively, of the original isospin-conserving sd 
interaction. We notice that in a fit to the Range III data, 
the charge-asymmetric part of the interaction increases 
up to 2.3-3.2% of the nuclear interaction. 

The Miller-Spencer parameterization and UCOM SRC 
schemes quench the Coulomb expectation values more 
than the AV-18 and CD-Bonn parameterizations (see Ta- 
ble II as an example). This is why the highest values 
of Xcoui in columns 2, 4, and 6 belong to UCOM SRC 
and Xcoui of the Miller-Spencer parameterization SRC are 
very close to them. At the same time, those parameteriza- 
tions result in the most negative values of Aq"^^ in columns 

2, 4, and 6 to compensate for the Coulomb efi'ect. 

b. Vcoui o,nd Vp combination. For the combina- 
tion of the Coulomb and Yukawa p-exchange type poten- 
tials as the isospin-symmetry breaking forces it should be 
noted that typical expectation values of Vp are about 3 to 
4 orders of magnitude smaller than the expectation val- 
ues of Vcoui- Therefore, small variations in the Coulomb 
strength (of the order of 1-2%) require a factor of up to 

20 variation in the corresponding strength Ap'^^ (e.g., the 
magnitudes of isovector p-exchange strengths range from 
4.4 to 100 for the USD interaction in the Range I data 
selection). So, the p-exchange potential strength is very 
sensitive to the SRC procedure. It is not surprising that 
the lowest absolute x'^p^ in columns 3, 5, and 7, corre- 
sponds to evaluations without taking SRC into account; 
however, with the presence of SRC, the lowest absolute 
X^p^ are from the CD-Bonn parameterization. The lowest 

Xcoui and, therefore, the highest absolute X^ff' in columns 

3, 5, and 7 belong to the Miller-Spencer parameterization 



SRC. The Xcoui values being closest to unity are from the 
CD-Bonn and UCOM SRC schemes for the three ranges 
of data (see also Fig. 3). 

Again we notice an increase of the isovector parameter 
Ap^'' in a fit to the Range III data. It was mentioned in 
Ref. [21], at that time it was not possible to conclude in 
that work whether the small asymmetry of the effective 
interaction is due to the original asymmetry of a bare NN 
force, or whether it was a radial wave function effect [69]. 
Our present results show that including more and more 
excited states lead to a bigger asymmetry value. Most 
probably, it is due to the radial wave functions effect, 
since asymmetry in the proton and neutron wave func- 
tions becomes larger in higher excited states. 

Regarding the ISPE's, their values stay consistent 
within certain intervals. Amazingly, the value of ^Q^g^j 
stays almost constant, without showing any dependence 
on the particular SRC approach, most probably, because 
it is the orbital which is most constrained by the data. At 
the same time, the value of £0)^3/2 depends somewhat on 

the SRC procedure. The highest value of £q|;3/2 column 
2 always corresponds to the Millcr-Spencer parameteriza- 
tion SRC. The values of £g)j3/2 ^-nd £^1^1/2 change much 
less for the Vcoui and Vp combination, than for the Vcoui 
and Vq combination (with the exception of the USDB in- 
teraction in the Range II data fit). As a general trend 
we notice a reduction of the values of ISPE's when we 
increase the number of data points in a fit. 

The values of parameters given in Table V lie outside 
the intervals obtained by Ormand and Brown who con- 
sidered the Vcoui and Vq combination. In particular, we 
get systematically lower values of the isotensor strength 
parameter Aq^"* , as well as lower values of eg|j5/2 ' oven for 
the Range I of data. 

The inclusion of nuclei from the middle of the sd shell, 
combined with the latest experimental data and with the 
newly developed approaches to SRC allowed us to con- 
struct a set of high-precision isospin-violating Hamilto- 
nians in the full sd shell-model space. They reproduce 
the experimental b and c coefficients with very low rms 
deviations, c.f. Table V. The ratios of the rms deviations 
of various SRCs, with USD, USDA, and USDB to the av- 
erage \b\ coefficients in sd-shell space are less than ~0.01. 
A few applications of these Hamiltonians are considered 
in the next sections. 



IV. THEORETICALLY FITTED h AND c 
COEFFICIENTS 

Theoretical b and c coefficients discussed in this sec- 
tion are obtained in a fit of the parameters of the charge- 
dependent Hamiltonian, consisting of the Coulomb inter- 
action (with UCOM type of the SRC) and Vq, on top 
of the USD interaction, to the experimental data from 
Range I. The obtained numerical values, as well as cor- 
responding experimental data, are given in Tables XII- 
XV for doublets, triplets, quartets, and quintetes, respec- 
tively. 

Before we begin the discussion, let us consider pre- 
dictions for b and c coefficients given by the uniformly- 
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FIG. 4. (Color online) Experimental and theoretical |6| coefficients in sd-shell nuclei plotted as a function of A. Minor x-axis 
tics are J states, which are arranged in an increment of 0.05 for every i step. 



charged sphere model [73]. In this approach, the total 
Coulomb energy of a nucleus is considered as a uniformly 
charged sphere of radius R = vqA^^^, 



E,, 



3e2 



5R 



Z{Z-l) 



(1 - A)T, + 



, (24) 



giving rise to the following expressions for a, b and c 
coefficients of the IMME [10, 74]: 

3e2 A{A - 2) 



20ro A3 
3e2 {A - 1) 



(25) 



5ro A3 
3e2 1 
or-o As 

where = 1.44 McV4m. Assuming Z{Z - I) k, Z"^ in 
Eq.(24), one can get an even simpler form of b coefhcients, 
namely. 



5ro 



(26) 



We will use these estimations and, in particular, the char- 
acteristic dependence on A in our analysis of the values 
obtained by a shell-model fit. 



A. Fitted h Coefficients 

Theoretical h coefficients are shown in Figs. 4 and 5 
in comparison with experimentally deduced values [62]. 
Overall, the deviations of b coefficients obtained in a shell- 
model fit from the experimental ones are less at the top 
and the bottom of sd-shell space than the deviations in 
the middle shell. The only exception is the \^ doublet of 
A ~ 39 for which the difference between theoretical and 
experimental b coefficients comes out to be 107.6 keV. 
However, if we refit the Hamiltonian parameters accord- 
ing to the smaller data range selected in Ref. [21] (the 
bottom and the top of the sd shell), this deviation for 

the doublet of A = 39 reduces to 49.8 keV. Thus the 
reason for a noticeable discrepancy for that point in a full 
sd shell Range I fit may be due to the inclusion of data 
from the middle shell. On the other hand, if we refit the 
parameters using extended data sets, Range //and Range 
III, this deviation reduces to 82 keV and to 0.7 keV, re- 
spectively. It is because the addition of more data points 
renormalizes the discrepancies of the fit. Although the 
inclusion of the b coefficient of the \^ doublet of A = 39 
reduces the quality of the fit, we retain it in the data set 
to adjust the ISPEs ef^ in Eq.(7). 

Let us remark that the quality of the fit is already 
somewhat pre-dctermined by the quality of the original 
isospin-conserving two-body interaction. For example, a 
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very accurate description of low-lying states in A = 35 
nuclei by the USD interaction leads to the values of the- 
oretical b coefficients of the A = 35 doublets which are 
close to the experimental ones (see Table XII). Another 
factor, the major factor, that influences the values of the 
obtained deviations is a characteristic property of the 
error- weighted least-squares fit. Experimental b coeffi- 
cients with very low error bars are favored in the shell- 
model fit and the corresponding theoretical b coefficients 
have typically rather low deviations. This is the reason 
why most of the lowest-lying multiplets' b coefficients are 
very close to experimental values. For example, the devi- 
ation between theoretical and experimental b coefficients 
of the mass A = 32 quintet, the best known quintet in the 
sd shell, is the lowest among the five quintets. Therefore, 
advances in mass measurements and nuclear excitation 
energies providing data points with low error bars may 
influence the data, which are dominant in adjusting the 
strengths of charge-dependent forces in the INC Hamil- 
tonian; in particular, data from the top and from the 
bottom of sd-shell space, which are used to calibrate the 
ISPEs. Similar magnitudes of deviations are obtained for 
other combinations of charge-dependent forces. For the 
USDA and USDB interactions (with either Vcoui + Vb or 
Vcoui + Vp, and with different SRC schemes), the devi- 
ations are a few keV higher than those obtained in the 
calculations with the USD interaction. 



As suggested by Eq.(26), we plot experimental and the- 
oretical b coefficients (obtained from a shell model fit) as 
a function of A^^^ in Fig. 5. It is evident that theoretical 
values are in remarkable agreement with the experimen- 
tal data. For comparison, we show b coefficients from the 
uniformly charged sphere model as well. Predictions of 
the latter reproduce well the trend of the b coefficients, 
however, they are about 500 keV off the experimental 
values when given by Eq.(25) and there is even a larger 
discrepency (about 800 keV) for a simplified form given 
by Eq.(26). Clearly, the ratio 1/Z for so?-shell space nuclei 
is not negligible with respect to Z. 

Although it is not seen in the scale of Fig. 5, in fact, 
the b coefficients in the sd shell-model space exhibit some 
regular oscillations around a straight line, representing 
the best fit of a linear function of A^^^ to the data points. 
This effect is discussed in Section IV C. 



B. Fitted c Coefficients 

The c coefficients obtained in the shell-model fit (see 
Tables XIII - XV) are ploted in Figs. 6 and 7. The 
discrepency between theoretical and experimental values 
for nuclei from the top and from the bottom of sd shell- 
model space are larger than for nuclei from the middle of 
the shell. Possible reasons for this have been mentioned 
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FIG. 6. (Color online) Experimental and theoretical c coefficients plotted as a function of A. Minor x-axis tics are J states, 
which are arranged in an increment of 0.05 for every i step. 



above, namely, it can be due to the larger experimen- 
tal error bars and/or lower accuracy of the correspond- 
ing isospin-conserving Hamiltonian to describe the energy 
levels. In Fig. 7, we plot the c coefficients as a function 
of A"^/'' as suggested by Eq.(25). Let us remark a few 
interesting features: 

• One easily notices a well pronounced oscillatory 
trend in the lowest-lying triplets' c coefficients con- 
nected by the solid line in Fig. 7. These values are 
always the highest or the lowest c coefficients, ex- 
cept for A = 20, 24, 28, and 32. This trend is also 
inherent to the corresponding experimental c coef- 
ficients [62]. 

• The first higher-lying triplets' c coefficients also ex- 
hibit regular oscillations, but of a smaller amplitude 
than those described above. The corresponding 
shell-model data points are connected by a double- 
dot-dashed line. 

• The other higher-lying triplets' and quartets' c coef- 
ficients lie somewhere in the middle part of the plot 
between maxima and minima of low-lying triplets' c 
coefficients without any particular behaviour. The 
quartets' c coefficients do not display any stagger- 
ing effect. The quintets' c coefficients connected by 
the dot-dashed line follow well the prediction of the 
uniformly charged sphere, c — ^A^^^'^ (dashed 



line). 

The theoretical c coefficients are seen to be in very good 
agreement with the experimental data. The uniformly 
charged sphere model describes well the overall trend of c 
coefficients, but it cannot predict the oscillatory behavior 
of the c coefficients. 



C. Staggering Behavior of b and c Coefficients 

The oscillatory effects in IMME b and c coefficients 
were noticed by Janecke in the 1960's, c.f Refs. [16, 
75], although at that moment the available experimental 
data was limited to T < 1 multiplets. Since then, a 
few analytical models have been proposed to explain the 
oscillatory effect. One of the approaches, proposed by 
Hecht [76], was based on Wigner's supermultiplet scheme. 
Another explanation was given by Janecke [16, 75] in the 
framework of a schematic approach to Coulomb pairing 
effects. 

In this section we revisit the staggering effect of the 
b and c coefficients of sd-shell nuclei based on a much 
more extended set of experimental data, which fully 
covers the lowest-lying doublets, triplets, quartets, and 
quintets, and we explore it theoretically using the con- 
structed empirical INC shell-model Hamiltonian. For the 
first time, we identify contributions of various isospin- 
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symmetry breaking terms to b coefficient (isovector en- 
ergy) and c coefficient (isotensor energy). 



1. Perspective of Empirical INC Hamiltonians 

To evidence a staggering pfienomenon, we plot the b 
coefficients obtained from experiment and from a shell- 
model fit for the lowest-lying doublets and quartets in 
sd-shell nuclei in Fig. 8(a). The oscillatory behavior of 
the b coefficients of doublets and quartets is clearly seen 
now. The data points form two families for A ~ 4n + 
1 and A = An + 3 multiplets lying slightly above and 
under the middle straight line, respectively. There is no 
staggering effect in the b coefficients of T = 1 triplets, 
c.f. Fig. 8 (d). This general behavior of the b coefficients 
of doublets, quartets and triplets agree with what had 
been noticed by Janecke [75] and by Hecht [76]. The 
quintets' b coefficients are known only for the lowest A = 
An multiplets and therefore we cannot discuss them on 
the same footing due to missing data. 

To magnify the effect of oscillations, we show devi- 
ations of the experimental and theoretical values from 
fitted middle lines (solid line for T = 1/2 multiplets 
and double-dot-dashed line for T = 3/2 multiplets) in 
Fig. 8(b) and Fig. 8(c). Interestingly, the oscillations 
of doublet b coefficients are of a higher amplitude com- 
pared to those of quartet b coefficients and they are in 



opposite direction. This tendency is naturally manifested 
in Wigner's supcrmultiplet theory [16, 77, 78]. As seen 
from these figures, the b coefficients obtained in a shell- 
model fit for doublets and quartets follow the experimen- 
tal trend extremely accurately, reproducing very precisely 
the general trend and the staggering amplitude. 

Since, the charge-dependent term in the INC Hamil- 
tonian is given by a combination of three components, 
KoidVcouh AqVo and ISPEs, J^i'^ty can explore what 
contribution from each component to the total b value 
is. The results are shown in Fig. 9 and Fig. 10 for the 
doublets' and the quartets' b coefficients, respectively. 

Qualitative analysis leads to rather similar conclusions 
for both doublets' and quartets' b coefficients. The isovec- 
tor Coulomb component is the main contribution to the 
staggering effects of the b coefficients of doublets and 
quartets (Figs. 9 and 10, respectively). It is interesting 
to note that the isovector charge-dependent term of nu- 
clear origin, \oV^'^\ produces the same oscillatory trend 
as that from the Coulomb force, but of a much smaller 
amplitude. The one-body contribution, however, does 
not produce any oscillations. This could be expected, 
since the staggering effect is due to the manifestation of 
the Coulomb contribution to the pairing. 

In general, the values of higher-lying multiplets' b co- 
efficients follow more and more smooth trends and the 
staggering gradually disappears. 

The general features of staggering have already been 
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FIG. 8. (Color online) Staggering efTect of the 6 coefficients of the ground-state doublets and of the lowest-lying quartets in 
sd-shell nuclei. Plot (a): |&| values of the ground-state doublets (left y axis) and the lowest-lying quartets (right j/-axis). For 
T = 1/2, the solid (blue) fine b = 0.7447A^''^ - 1.2551 (MeV), the dashed (black) line 0.7442A^''^ - 1.1987, and the dotted-dash 
(black) fine b = 0.7463^^/^ - 1.3329 (MeV), represent best fits to experimental values for the all A, A = 4n -I- 3 and A = 4n + 1 
isobars, respectively. For T = 3/2, the double-dot-dashed fine 6 = 0.7441A^/^ - 1.2547 (MeV) is obtained by a fit to all 
experimental b coefficients. Plots (b) - (d): deviations of b coefficients from the respective fitted to all data points lines. The 
experimental 6 values are represented by (blue) triangles and b obtained in a shell-model fit are shown by (red) squares. 



discussed in section IV B. In Fig. 11 we plot separately 
the c coefficients of the lowest-lying triplets (the upper 
part of the figure) and the lowest-lying quartets and quin- 
tets (the lower part of the figure) in sd-shell nuclei. The 
c coefficients obtained in the shell-model fit reproduce 
the experimental values very precisely (with the largest 
deviation of about 15 keV, see also Tables XIII - XV). 

The experimental and shell model fitted c coefficients 
of triplets clearly form two distinct families of multiplcts 
for A = An and A = An + 2 nuclei, respectively. However, 
no oscillations can be noticed for T = 3/2 multiplets. 
The c coefficients of quintets are known only for A = An 
multiplets which follow a quite smooth trend with mass 
number. 

Contributions of different terms of the charge- 
dependent Hamiltonian to the lowest-lying triplets' c co- 
efficients are shown in Fig. 12. One can see that the 

(2) 

isotensor Coulomb force V^^^i plays the major role. Fur- 

(2) 

thermore, the plot also indicates that V^o„/ alone does 
not reproduce the magnitude of the experimental c coef- 
ficients. For the A = An + 2 family, the deviation is about 
~40 keV, while for the A = An family, it is around ~5 
keV. This indicates that the Coulomb interaction should 
be supplemented by another two-body interaction of nu- 



clear origin, which we model as Vq' (or Vp ) in this 
paper and which perfectly fulfills its task. The contri- 
bution of the empirical isotensor nuclear interaction re- 
sults in the same oscillatory trend as that of the isotensor 
Coulomb component, with the values being of about ~40 
keV for A ~ An multiplets and ~5 keV for A = An + 2 
multiplets (with a negative value for the lowest ^ = 36 
triplet). Thus, the experimental values of c coefficients 
are perfectly reproduced. 

A similar decomposition of the theoretical c coefficients 
for quintets is given in Fig. 13. As has been already men- 
tioned, the data on A = 4n -I- 2 multiplets is required 
in order to establish the existence of the staggering ef- 
fect. It is seen, however, that the contribution from the 
isotensor nuclear force to quintets' c coefficients shows 
some noticeable oscillatory effect between A = 8n and 
A = 8n + A multiplets. It may possibly change to the 
staggering characteristics for triplets' c coefficients when 
data on A ^ An + 2 becomes available. 

The c coefficients of high-lying multiplets are system- 
atically known only for triplets. As it was mentioned in 
the previous section, the first high-lying triplets' c coef- 
ficients oscillate with a smaller amplitude, while c coeffi- 
cients of other high-lying multiplets follow a more or less 
smooth trend. This is probably related to the destroying 
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FIG. 9. (Color online) Contributions of the various charge- 
dependent forces to doublet (T = 1/2) b coefficients. The |6| 
values are plotted as a function of A^^^ . The total |6| values 
refer to the left y-axis, while contributions from V^^li, from 

refer to the right 



Vr 



and from the summed ISPE e. 



(1) 



y-axis. The contribution from V^'^' is multiplied by 20 



of the pairing effects with increasing excitation energy in 
nuclear systems. 

Very similar trends and exactly the same conclusions 
can be inferred if other sd-shcU model interactions are 
used instead of USD, or other charge-dependent Hamil- 
tonians (with other SRC schemes). This proves the ro- 
bustness of the effects described above. 



2. Jdnecke's Schematic Model 

Janecke's model [16, 75] is based on an approximate 
formula for the Coulomb energy of valence proton(s) 
outside a closed shell, which was proposed by Carlson 
and Talmi [79]. In order to match the trend of the to- 
tal Coulomb energy of a nucleus as represented by the 
IMME, Janecke replaced the Coulomb pairing term [79] 
by a quadratic term in [75]. As a result, one can de- 
duce the following expressions for isovector, E^^l^^, and 
isotensor, E^J^^, contributions. 



E 



(1) 



Ie,a 



E2 + ^e^ . 



(27) 



19 21 23 25 27 29 31 33 35 37 39 
A 

FIG. 10. (Color online) Contributions of the various charge- 
dependent forces to the lowest-lying quartet (T = 3/2) & co- 
efficients. Refer to Fig. 9 for further description. 



and 



2uE^) 



(28) 



where the energies Ei are related to one and two-body 
electromagnetic interactions, while /i and v are some pa- 
rameters. In Ref. [75], assuming an independent-particle 
model with four-fold degenerate orbitals, Janecke could 
estimate a probability for the number of proton pairs to 
occupy the same orbital and, thus, he could deduce a 
contribution to the Coulomb energy for a given A and T. 
He obtained the following parameterization for /i and v 
values: 



M = 



1 

2 

5'- 



(-i)-v 



2-T 



2T 



A- even 



A-odd 



(29) 



and 



^ 1 



1 

4T 
1 

4T 



2T-1 



A- 



A-odd, T > i . 



(30) 



As was remarked in Rcf. [75], the coefficients Ei with 
i = 1,2, and 3, are related to the expectation value of 1 /r. 
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FIG. 11. (Color online) Staggering effects of c coefScients 
of the lowest-lying triplets (upper panel), and quartets and 
quintets (lower panel). See text for discussion. In the lower 
panel, the dot-dashed (black) line, c = 651.6A~^''^ + 2.6 (keV), 
is an unweighted fit to the experimental c coefficients. The 
solid (blue) line, c = 520. 7^-^/^-^44.1 (keV) , is an unweighted 
fit to the c coefficients obtained in a shell-model fit. 



because the average distance between protons should in- 
crease with the nuclear volume. Therefore, if we assume 
that Ei = Ei/A^/"^, with Ei being constant values, differ- 
ent from one shell to another shell, the isovector Coulomb 
energy will become a linear function of A^^^ and the 
isotensor energy will be a linear function of A^^^^ (since 
El is the leading term in expressions (Eq.(27)-Eq.(28)) . 

a. Isovector Coulomb Energies. Using Eq.(27) and 
the respective fj, values, we may derive the isovector 
Coulomb energies for doublets, triplets, quartets and 
quintets as 



E 



(1) 

coul (T=l/2) 
(1) 



lEiA + E,+ '-E, + {~ir-^^y'^, 

(31) 



E 



E 
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(1) 



^EiA + E2 + ^E3 



(32) 



lEiA + E,+ '-Es + (-ir~'^^'^, 

(33) 



E 



coul (T=2) 



^EiA + E2 + ^E3 



(34) 



FIG. 12. (Color online) Contributions of the various charge- 
dependent forces to the lowest-lying triplet c coefficients. All 
c values refer to the left j/-axis, while the contributions from 



V}1, and V, 



(2) 



refer to the right y-axis. 



respectively. From the last terms of Eq.(31) and Eq.(33), 
determining the amplitude of the oscillations of the b 
coefficients of doublets and quartets, we see that the 
amplitude of quartets' c coefficients is predicted to be 
three times smaller than that for doublets. Eqs. (32) and 
(34) indicate that no oscillatory behavior is expected for 
triplets' and quintets' isovector Coulomb energies (or b 
coefficients). 

b. Isotensor Coulomb Energies. From Eqs. (28) and 
(30), we can obtain the isotensor Coulomb energies for 
triplets, quartets and quintets as 



1 



1 



1 - (-l)^/2i?3 



(2) 1/1 

^coul{T=S/2) = 6 V ^ ^ 3^^ 



^coul(T=2) ~ 6 I 4 



1 - ^(-i)(^-^'/^£;3 



(35) 
(36) 

)• 

(37) 



respectively. The last term of Eq.(35) and Eq.(37) shows 
that triplets' and quintets' c coefficients exhibit regu- 
lar oscillations as a function of A, with the amplitude 
for triplets being six times larger than that for quintets 
(c = '^^coui'' ^"^^ Eq.(9)). Eq.(36) shows that the quar- 
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FIG. 13. (Color online) Contributions of the various charge- 
dependent forces to the lowest-lying quintet c coefficients. All 
c values and the contribution from V^^^j refer to the left j/-axis, 
while the contribution from Vq'^' refer to the right y-axis. 

(2) 

tets' isotensor Coulomb energy -E'com/ (t=3/2) predicted 
to be a constant which may vary from one shell to an- 
other. Hence, an oscillatory behavior is not predicted for 
quartets' c coefficients. 

Performing a linear fit to the experimental h coefficients 
for the lowest-lying doublets, we have determined the val- 
ues of El = 487 kcV, E2 = 1199 keV, and = 134 keV 
for sc?-shell nuclei. The value of E3 deduced from the fit 
to b coefficients predicts the ^i?3 = 68 keV amplitude for 
c coefficients in T = 1, which is in very good agreement 
with the experimental value. Analysis of staggering in 
other model spaces and the values of Ei coefficients will 
be published elsewhere [78]. 



V. MASSES AND EXTENSION OF THE IMME 
BEYOND THE QUADRATIC FORM. EXAMPLE 
OF THE A = 32 QUINTET. 

The fit and the analysis in Section IV C are based on 
the assumption of a quadratic form of the IMME, which 
is a very good approximation, valid at present for the 
majority of experimentally measured isobaric multiplcts. 
However, some experimental cases evidence the breaking 
of the quadratic IMME (e.g., see Refs. [62, 64] and refer- 



ences therein, as well as Refs. [80-83]). We consider here 
an extended IMME up to a quartic form. 



M{a, T, T,) =a(a, T) + 5(a, T)T, + c(a, T)T^^ 

+d(a,T)T3 + e(a,T)^,^ (38) 



with possible non-zero d and/or e coefficients. These 
higher-order terms in can be due to the presence of 
isospin-symmetry breaking three- (or four-body) inter- 
actions among the nucleons [84], and/or may arise due 
to the isospin mixing in excited states of isobaric multi- 
plcts with nearby statc(s) of the same J'^, but different 
T value. Early theoretical estimations for quartets pre- 
dicted typical d coefficients to be of the order of wl keV 
[85-87] (see also discussion in Ref. [74]). To probe such 
low values, recent experimental advances become crucial 
in providing precise mass measurements of quartets and 
quintets. At present, relative mass uncertainties as low 
as 10-* - 10-9 are reached, see e.g. Refs. [80, 82, 83, 88]. 

In the shell model the direct evaluation of absolute 
binding energies is possible with the isospin-conscrving 
Hamiltonian, provided that a certain algorithm is fol- 
lowed in the subtraction of empirical Coulomb energies 
from experimental binding energies used in the fit. Then, 
the subtracted Coulomb energy should simply be added 
to the shell-model binding energy to get the full theoreti- 
cal binding energy of a nucleus. In fitting the USD inter- 
action, the subtraction of the Coulomb energy has been 
done in a kind of average way [89, 90]. In particular, an 
unknown amount of residual isoscalar Coulomb energy 
may remain in the charge independent nuclear Hamilto- 
nian [89, 90]. Adding an INC term in the Hamiltonian 
requires the precise knowledge of the isoscalar Coulomb 
contribution and this prohibits the evaluation of abso- 
lute binding energies [28]. In spite of this fact, we can 
still well describe theoretical mass differences of isobaric 
multiplcts, which is sufficient to study the 6, c, d, and e 
coefficients of the IMME. The a coefficient, however, re- 
mains undetermined. To explore theoretically the validity 
of the quadratic, cubic or quartic forms of them IMME 
in a given quintet, we use the results of the exact diag- 
onalization of the INC Hamiltonian, i?iNCj constructed 
in the present work. In this way we obtain theoretical 
mass differences for a given isobaric multiplct and then 
we fit them with a quadratic, cubic or quartic form of the 
IMME to find the best 6, c, d, and/or e coefficients. 

As an example, here we consider in detail the lowest 
0+ quintet in A = 32. 

Various experimental determinations of the lowest T = 
2 masses in ^ = 32 [80, 82, 83] point towards the presence 
of a non-zero d coefficient in the IMME (see Table VIII 
later in this section). 

Using Eq.(38), we can express the IMME a, 6, c, d, and 
e coefficients in terms of the mass excesses of a given 
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TABLE VI. Mass differences and mass summations of A/_2 and M2; and Mi, M_i, and Mq. 





M_2 - A/2 

(kcV) 


Ml - M_i 
(keV) 


Afi + A/_i - 2Afo 
(koV) 


A/2 + A/_2 - 2A/o 
(keV) 


Exp. values 
quoted in Ref. [91] 


-21877.48 


-10944.24 


414.22 


1657.26 


Theoretical values 
from Ref. [91]: 
USD 


-21669.83 


-10837.25 


418.11 


1673.09 


USDA 


-21669.62 


-10836.63 


404.98 


1653.43 


USDS 


-24802.19 


-12402.79 


417.24 


1667.38 



-10943.63 413.91 1657.05 

-10927.87 414.87 1660.55 

-10927.30 404.55 1649.76 

-10927.93 416.24 1664.32 



Deduced from Table XV -21877.29 
Present work^: 

USD -21857.96 

USDA -21858.35 

USDS -27858.34 



Present calculations use Vcoul (UCOM) and Vq combination. 

T = 2 quintet, 
a = Mo , (39a) 
b = ^ [(A/_2 - M2) + 8(Afi - A/_i)] , (39b) 

c = ^ [16(A/i + A/_i) - (A/2 + A/_2) - 30A/o] , (39c) 

d=^ [(A/2 - A/_2) + 2(A/i - A/_i)] , (39d) 

e = ^ [-4(A/i + A/_i) + (A/2 + A/_2) + 6A/0] . (39e) 

Here, we have shortened the notation for a, &, c, d, and 
e coefficients and the notation for mass excess of each 
member (Mt^^, = Mi, i = -2,-1,0,1,2). Eqs. (39b) 
and (39d) show that h and d coefficients are related to 
the differences (A/_2 - A/2) and (A/i - A/_i). Note that 
h and d are not linked to a. Meanwhile, c and e are de- 
fined by the a-subtracted sums of A/2 and A^_2 and the 
sum of Ml and M-i. These coefficients are also inde- 
pendent of a (then a enters in each mass member and 
cancels in the expressions Eq.(39c) and Eq.(39e)). This 
set of relations is kept for 4-parameter least-squares fits 
to the cubic IMME (Eq.(38) with e = 0), or to the quar- 
tic IMME (Eq.(38) with d = 0), or is solved exactly in 
the case of the full quartic IMME (both d and e arc non- 
zero). In our theoretical analysis, we assume that every 
input mass excess has the same uncertainty, e.g., ±1 keV. 

Table VI summarizes mass differences (or sums) of ibTj 
multiplet members as obtained from the experimental or 
theoretical mass excesses. We have performed calcula- 
tions using all the USD, USDA and USDB interactions 
and the combination of Vcoui (with UCOM) and Vq as an 
INC term with the parameters found by the fit (c.f. sec- 
tions II and IIIC). The obtained results (the lower part 
of Table VI) are compared with the recent analysis of Sig- 
noracci and Brown [91], who performed a similar study. 



but using the INC Hamiltonian parameterization from 
Ref. [21] (the upper part of the same table). It is seen 
that the mass differences (A/_2 - A/2) and (A/i - A/_i) 
obtained in the present work are systematically closer to 
the experimental values than those of Ref. [91]. These 
are exactly the key figures which determine h and d coef- 
ficients. 

Table VII shows theoretical IMME 6, c, </, and e co- 
efficients obtained for each set of mass differences by a 
least-squares fitting procedure assuming all uncertainties 
of the theoretical mass excesses of A = 32 to be 1 keV. 
The present results (the lower part of the table) are com- 
pared with the results of Signoracci and Brown (the upper 
part of the table). Two slightly different sets of experi- 
mental mass excesses are taken from Ref. [91] (the first 
entries in the upper and lower parts of Table VII). 

As seen from Eq.(39b) and Eq.(39b), the presence of 
the d coefficient adjusts the respective h coefficient in the 
fit. Theoretical b coefficients in the third and the fifth 
column are the same, since the d coefficient is not con- 
sidered in the corresponding fits. Similarly, b coefficients 
in the fourth column and the last column are the same, 
because the d coefficient is included in those fits. A sim- 
ilar situation holds for the c and e coefficients, which are 
determined by the a-removed sum of the mass excesses of 
= ±1 and Tz = ±2 isobaric members of the multiplet 
as follows from Eq.(39c) and Eq.(39e). 

Before we discuss any evidence for non-zero c? or e co- 
efficients, let us compare the values of the corresponding 
b and c coefficients. As seen from Table VII, b coefficients 
obtained in the present work reproduce much better the 
experimental values compared to the results of Ref. [91]. 
In particular, we get all deviations smaller than 10 keV, 
while the calculations of Ref. [91] result in much larger de- 
viations of -50 keV (USD, USDA) to -730 keV (USDB). 
This is due to the fact that the corresponding mass dif- 
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TABLE VII. Comparison 


of b, c, d, and e 


coefficients of the A = 


32, J'^ = 0+,r = 2 quintet. 








b, c 


b, c, d 


b, c, e 


b, c, d, e 






(keV) 


(keV) 


(keV) 


(keV) 


Exp. values 


L 




-5471.9 (3) 


-5473.1 (3) 


-5471.1 (3) 


-5473.0 (5) 


quoted in Ref. [91] 


c 


208.6 (2) 


207.2 (3) 


205.5 (5) 


207.1 (6) 




d 





0.93 (12) 





0.92 (19) 




e 








0.61 (10) 


0.02 (16) 




2 / 

X /n 


32.14 


0.01015 


22.63 




Theoretical values 












in Ref. [91]: 












USD 


b 


-5417.7 


-5419.0 


-5417.7 


-5419.0 




c 


209.1 


209.1 


209.0 


209.0 




d 




0.39 




0.39 




e 







0.03 


0.03 






1.089 


0.005491 


2.172 


— 


USDA 


b 


-5417.6 


-5418.6 


-5417.7 


-5418.6 




c 


207.3 


207.3 


201.1 


201.1 




d 




0.30 





0.30 




e 








1.40 


1.40 




X /n 


8.680 


16.04 


1.318 


— 


USDB 


b 


-6200.7 


-6201.7 


-6200.7 


-6201.8 




c 


208.4 


208.4 


208.7 


208.7 




d 




0.28 




0.28 










-0.07 


-0.07 




xVn 


0.5931 


0.03703 


1.1492 





Present work^: 

Exp. values 

taken from Table XV 



xVn 



-5471.85 (27) 
208.55 (14) 



32.15 



-5472.83 (29) 
207.12 (23) 
0.89 (11) 

0.1035 



-5470.45 (29) 
204.92 (23) 

0.69 (11) 
13.80 



-5472.64 (68) 
206.89 (75) 
0.83 (22) 
0.06 (19) 



USD 



b 
c 
d 

xVn 



-5464.38 
207.59 



0.2563 



-5463.75 
207.59 
-0.19 

0.01685 



-5464.38 
207.39 

0.045 
0.4958 



-5463.75 
207.39 
-0.19 
0.045 



USDA 



b 
c 
d 



xVn 



-5464.40 
206.78 



7.805 



-5463.34 
206.78 
-0.31 

14.21 



-5464.34 
200.96 

1.31 
1.400 



-5463.34 
200.96 
-0.31 
1.31 



USDB 



b 
c 
d 

xVn 



-5464.13 
208.03 



0.03538 



-5463.90 
208.03 
-0.07 

0.006007 



-5464.13 
208.15 

-0.03 
0.06475 



-5463.90 
208.15 
-0.07 
-0.03 



^ Present calculations use Vcoul (UCOM) and Vq combination. All jn are given in 4 significant figures. 



ferences (M_2 — M2) deviate from the experimental value 
by about 207 keV (USD, USDA) and 2924 keV (USDB), 
while the mass differences (Mi — M„i) are different from 
the experimental value by about 107 keV (USD, USDA), 
and 1458 keV (USDB) (see Table VI). The present INC 
Hamiltonian produces mass differences which deviate at 
most by 20 kcV from the experimental values and thus 
results in very close h values. This discrepancy should be 



kept in mind when comparing the values of the predicted 
d coefficient with the experimental value. 

Let us remark that the theoretical h coefficients calcu- 
lated via an exact diagonalization almost coincide with 
the h coefficients listed in Table XV, which were obtained 
in a fit (within perturbation theory). That means, the 
perturbation theory used in section II provides a very 
good approximation to the h coefficients. 
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FIG. 14. (Color online) Deviations of experimental or the- 
oretical masses of the A = 32 quintet from the correspond- 
ing quadratic IMME fit. (Purple) squares are quoted from 
Ref. [80]. (Purple) filled squares are quoted from Ref. [82]. 
(Green) circles are quoted from Ref. [83], Table VIII, set A; 
(green) filled circles are from set B; up (blue) triangles are 
from set C; up (blue) filled triangles are from set D; down 
(black) triangles are from set E; down (black) filled trian- 
gles are from set F. The recent theoretical work of Signoracci 
and Brown [91] is presented as (red) pentagons; whereas (red) 
filled pentagons are present calculation. 




FIG. 15. (Color online) Deviations of experimental or theo- 
retical masses of the A = 32 quintet from the corresponding 
cubic IMME fit. Refer to Fig. 14 for description. 



Overall the c coefficients predicted by both models are 
close to experimental values, with a maxium ~2 keV for 
the present results and even ~1 kcV for the values of 
Ref. [91]. 

Each set of IMME coefficients in Table VII is ended by 
the x^/n value characterizing the quality of the fit. It is 
seen that all calculations agree well with the experimental 
conclusion that the cubic form of the IMME describes 
best the nuclear mass trend of the lowest 0+ quintet in 
A = 32, since it produces the lowest x^/n value (with 
the exception of the USDA interaction, see explanation 
below). The quartic IMME with c? = is worse than the 
cubic one (again, except for the prediction of the USDA 
interaction) . 

To illustrate this effect, we plot in Fig. 14 and Fig. 15 
deviations of nuclear mass excesses from the best IMME 



fit values, assuming a quadratic and a cubic form of the 
IMME, respectively. These figures include different ex- 
perimental data sets and two different theoretical calcu- 
lations of mass excesses (Ref. [91] and present work). It 
is obvious that the best fit is produced by a cubic form 
of the IMME (Fig. 15). 

The values of the corresponding d coefficient are, how- 
ever, different in experimental and theoretical analysis. 
The experimental value ranges from 0.51 keV to 1.00 keV 
for various sets of experimental data (see Table VIII). 
Taken the adopted values of experimental mass excesses, 
we get dexp = 0.89(11) keV (from our recent compila- 
tion [62]). At the same time, theoretical values obtained 
from the USD interacton are dth = 0.39 keV [91] and 
dth = -0.19 kcV (present result). The USDB interac- 
tion results in much smaller x^/*^ values for all fits, with 
the minimum again for a cubic form of the IMME. The 
corresponding d coefficients are dth = 0.28 keV [91] and 
dth = —0.07 keV (present result). 

Let us remark that although the d- values of Ref. [91] 
are closer to the experimental one, there is an essential 
discrepency in their theoretical b coefficients, especially 
for the USDB interaction. At the same time, although 
being in better agreement for b coefficients, our calcula- 
tions point towards a negative value of the d coefficient. 
Wc think that this is due to a peculiarity of the fit, since 
the sign of the d coefficient is determined by a ratio of 
mass differences (M_2 -M2)/(Mi -M_i) (see Eq.(39d)). 
To get a zero value for the d coefficient, this ratio should 
be equal to 2. The ratio we get with the USD interac- 
tion is 2.00019, resulting in a negative d value, while the 
USD calculation of Ref. [91] produces a ratio of 1.99957, 
which is closer to the experimental ratio of 1.99902 and 
both producing a positive d value. We believe that ex- 
perimental data of high precision can certainly help to 
refine theoretical predictions. 

As was shown by Signoracci and Brown [91], the shell- 
model value of the d coefficient relates to the degree of the 
isospin-mixing in the initial and final nuclear states. For 
this reason, the calculation with the USDA interaction 
should be taken with caution, since accidentally there is 
a very closely lying state 0+, T = 1 in the vicinity of the 
IAS 0+, T=:2 state in ^^ci [91]. We arrive at a similar 
conclusion while adding our parameterization of the INC 
Hamiltonian to the USDA interaction. 

Apparently, the effects discussed in this section require 
a very high precision of relevant experimental data and 
theoretical accuracy. In particular, the essential uncer- 
tainty of "^^Ar mass may affect both experimental con- 
clusions and theoretical description. Hence, a direct re- 
measurement of the ^^Ar mass would shed light on this 
issue. 



VI. ISOSPIN-SYMMETRY BREAKING AND 
FERMI BETA DECAY 

A. Superallowed Fermi beta decay 

The absolute value of a superallowed O"*" — )> 0+ 
^-decay can be deduced from the experimental ft value 
applying various corrections (see Ref. [5] and references 
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TABLE VIII. Comparison of theoretical d coefficients with experimental values for the A = 32 quintet. 



ExpGrimGntfl.! cLiid TliGorGticcil Works 


d (keV) 


V / 77 .-1 

\ I ".quaar. 


'i 1 

\^ j "'Cubic 


Triamhak S pf nl 

_L1 IdlllUCLJV O. Ltd. OW 




6.5 


0.77 


Kwiatkowski S. et. al. [82] 


1.00 (9) 


30.6 


0.48 


OGl ri iiUIll JXtiliKcLiliGll jt.. Hi. (IL. lOOl 


U.OZ l^-i-^J 


Q Q 


U.oD 


Set B Ibid. 


0.60 (f3) 


f2.3 


0.3f 


Set C Ihid. 


0.90 (f2) 


28.3 


0.002 


Set D 


f.OO (fS) 


30.8 


0.09 


Set E Ihid. 


0.5f (f5) 


6.5 


0.74 


Set F /birf. 


0.62 (f6) 


8.3 


0.28 


Signoracci A. & Brown B. A. [9f] (USD)'' 


0.39 


f.09 


0.005 


Present work^ 


-0.f9 


0.26 


0.02 



^ X'^/"quadr. ^nd X^/"cubic of Refs. [80, 82, 91] are calculated in present work, moreover X'^/"quadr. and X^/"cubic of theoretical works 
are calculated by assuming an uncertainty of ±1 keV for every mass excess. 



TABLE IX. Comparison of 5im values. 



Nuclear Parent 
Hamiltonian Nucleus 


Present Work^ 


Previous Work 


w/o SRC 


UCOM 


Jastrow-type SRC function 
Argonne V18 CD-Bonn Miller 

Spencer 


Ormand & 
Brown 


Towner & 
Hardy'' 


USD 






















0.021 


0.022 


0.021 


0.021 


0.023 


0.017^= 


0.010 (10) 




26m 


0.011 


0.012 


0.012 


0.011 


0.013 


0.01=^ 


0.025 (10) 




2«Si 


0.046 


0.046 


0.046 


0.046 


0.046 


0.028^^ 


0.022 (10) 




30 g 


0.026 


0.028 


0.027 


0.025 


0.030 


0.056^^ 


0.137 (20) 




34C1 


0.037 


0.037 


0.036 


0.036 


0.036 


0.06^^ 


0.091 (10) 




34 Ar 


0.005 


0.006 


0.006 


0.006 


0.007 


0.008" 


0.023 (10) 


USDA 




















^^Mg 


0.021 


0.022 


0.021 


0.020 


0.024 








26m 


0.012 


0.013 


0.013 


0.012 


0.015 








2«Si 


0.041 


0.042 


0.042 


0.041 


0.043 








30 g 


0.020 


0.022 


0.021 


0.020 


0.024 








34C1 


0.031 


0.032 


0.031 


0.031 


0.032 








34 Ar 


0.006 


0.007 


0.007 


0.006 


0.008 






USDB 




















^^Mg 


0.021 


0.022 


0.021 


0.021 


0.023 








26m 


0.012 


0.013 


0.013 


0.013 


0.015 








^«Si 


0.044 


0.045 


0.044 


0.044 


0.046 








30 g 


0.025 


0.026 


0.025 


0.024 


0.028 








34C1 


0.036 


0.036 


0.035 


0.035 


0.035 








34 Ar 


0.005 


0.005 


0.005 


0.005 


0.006 







Present calculations use strength parameters from Table V corresponding to the Vcoui plus Vq combination. 
Unsealed &c\ from Table III of Ref. [24]. 

The <5iM obtained in the present work using the INC Hamiltonian of Ref. [21] and without truncation. 
The (5iM from Table I of Ref. [23]. 



therein), i.e. 

.9^t= /t(l + 5{,,)(l + <5Ns-fc) 
K 

= G2.|A/^o|2(i + AX) ' 

where A)^, 5^, and (5ns are transition-independent, 
transition-dependent and structure-dependent parts of 
the radiative correction, respectively, while be is 
the nuclear structure correction due to the isospin- 
symmetry breaking in the parent and daughter nuclear 



states. Other constants in this expression are K ~ 
2iT^hln2/{mec'^f = (8120.278±0.0004) x lO'^^ GeV-^s, 
Gv is the vector couphng constant for a semileptonic 
weak process, while 

\Mfo\^ = \{iJf\T+\i;,)\' =T{T + l)-T,,T,f 

is the value of the Fermi matrix element squared in the 
isospin-symmetry limit, which in the case of T = 1 emit- 
ters is |AfFoP = 2. 

Provided the absolute value of superallowed 0+ — 
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TABLE X. Comparison of values. 



Parent 
Nucleus 


Exp. ft (s) 


Sim 


Corrections (%) 

Sro 5h, 


3ns 


Present 
work 


(s) 

Towner & 
Hardy^ 


22 Mg 


3052 (7) 


0.0216 (9) 


0.370 (20) 


1.466 (17) 


-0.225 (20) 


3077.6 (72) 


3077.6 (74) 


26 ^ym 


3036.9 (9) 


0.0120 (8) 


0.280 (15) 


1.478 (20) 


0.005 (20) 


3072.9 (13) 


3072.4 (14) 


3*C1 


3049.4 (12) 


0.0363 (5) 


0.550 (45) 


1.443 (32) 


-0.085 (15) 


3072.6 (21) 


3070.6 (21) 


34 Ar 


3053 (8) 


0.0060 (4) 


0.635 (55) 


1.412 (35) 


-0.180 (15) 


3070.7 (84) 


3069.6 (85) 



^ Values quoted from Table 4 in Ref. [-5]. 



3090 



s3 3080 
-a 3070 
^ 3060 
^ 3050 
3040 



3030 



3090 



-a 3080 
Cj 3070 



3060 



Exp, ft 



74 



Rb 



64 



Ga 
i 



22 



Mg 



34 



Ar 



34 



'-C "o 



a 



'26 



'Al" 



10 
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Present work 
Towner and Hardy (2010) 
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FIG. 16. (Color online) Comparison of values in sd-shell space. The upper figure depicts experimental ft values. The 
ft value of ^^K"^ is not shown, it is 3051.9 (10) s. The lower figure shows ^t values. Horizontal (grey) strip is 1 standard 
deviation according to Ref. [5]. 



0+ /3 transitions is constant as stated by the CVC hy- 
pothesis, the Gv value ean be deduced from Eq.(40) and 
then the Vud CKM matrix element can be obtained from 
comparison of Gv with the vector coupling constant ex- 
tracted from the muon decay. 

The isospin-symmetry breaking correction is defined as 
a deviation of the realistic Fermi matrix element from its 
model-independent value: 

\Mf\^ ^\Mfo\^{1~Sc). (41) 

Within the shell model, the initial and final nu- 
clear states represent a mixing of many-body spheri- 
cal harmonic-oscillator configurations. In practice, the 
isospin-symmetry breaking correction to \Mfq\'^ is usu- 
ally separated in two terms: Sc = Sim + ^ro (we adopt 
here the notations of Ormand and Brown [20, 23]; in the 



work of Towner and Hardy [5, 24, 92, 93], these terms 
are referred to as Sci and Sc2, respectively). The first 
term, 6im-, is a correction to the Fermi matrix element 
MpQ from the isospin-symmetry breaking in the config- 
uration mixing of the spherical harmonic-oscillator basis 
functions (isospin-mixing correction). This is obtained 
via the diagonalization of an effective INC Hamiltonian 
within the valence space. The second term, 5ro, arises 
in the calculation of transition matrix elements due to 
the non-unity of the radial overlap of proton and neu- 
tron wave functions (radial-overlap correction). To get it, 
one has to replace the harmonic-oscillator single-particle 
wave functions by more realistic spherically-symmetric 
wave functions obtained from a better suited finite-well 
plus Coulomb potential (to account for the isospin non- 
conservation outside the model space). 
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TABLE XI. Values of Sim for Fermi 13 transitions between non-analog states (from the lowest 0"'',r — 1 state of '^^Ar or '^^Cl to 
the first excited 0^, T = 1 state of the corresponding daughter nuclei). 



Nuclear Present Work^ Previous Work 

Hamiltonian 

w/o SRC UCOM Jastrow type SRC function 

Argonnc V18 CD-Bonn Miller Strengths from 
Spencer Ref. [21] 



3*Ar 

USD -37.97 -8.843 -29.69 -48.19 -27.44 -361.2 -3.5 

USDA -127.8 -101.2 -118.9 -140.2 -71.25 — — 

USDB -134.1 -120.8 -134.5 -142.9 -106.1 — — 



Ormand & 
Brown [20] 



Sim X 10" : 
^^Cl 

USD -3.103 -3.098 -3.061 -3.071 -3.061 -3.960 -2.7 

USDA -2.734 -2.730 -2.697 -2.712 -2.695 — — 

USDB -3.120 -3.115 -3.076 -3.088 -3.071 — — 



The 5iM values of present work arc given in 4 significant figures. 



In the present paper, we present calculations of the 
isospin-mixing corrections Sim to the experimental ft val- 
ues for 0+ 0+ /? transitions in sd-shell nuclei. Although 
this correction is known to be quite small (from ~0.01% 
to ~0.1%), we can still see noticeable changes to the abso- 
lute J^t values compared to the existing evaluation, based 
on the calculations of Towner and Hardy [5] . 

The values of i5im obtained from the USD, USDA and 
USDB interactions with Vcoui + Vq INC Hamiltonian with 
all possible approaches to the SRC in Coulomb TBME's 
are summarized in Table IX. In the last two columns, we 
give for comparison the values from the previous work by 
Ormand and Brown [23] and the most recent results of 
Towner and Hardy [24]. 

As seen from the table, in spite of different interac- 
tions used, our calculations lead to very consistent values 
for various isotopes. Only for ^^Si and '^'^S, the USDA 
interaction results in somewhat smaller values of correc- 
tions than those predicted by the USD or USDB. The 
dependence on one or another approach to the SRC is 
also marginal. For most isotopes, the Miller and Spencer 
parameterization produces slightly larger values of cor- 
rections. 

If we compare our results with those obtained by Or- 
mand and Brown, we see that for ^^Mg, 26m^j 34^^. 
there is a good agreement. However, clear differences can 
be noticed for ^^Si, '^°S and ^"'Cl. Wc get a larger value 
for the former and almost twice smaller values for the two 
latter cases. Regarding the calculation of Towner and 
Hardy (last column), only for the two lightest emitters 
(^^Mg and ^^™A1) our values are close to their range (un- 
certainties included). For the rest of the nuclei, there is 
an essential discrepancy — we get a larger value for ^^Si 
and much smaller values for the heavier emitters. The 
main difference in the approaches is that the authors of 
Ref. [24] adjusted INC Hamiltonian strength parameters 
separately for each considered multiplet (case by case) to 
reproduce the isobaric mass splitting, while in our work 



and that of Ref. [23] a global parameterization for sd-shell 
nuclei has been exploited. 

In Table X we apply the calculated values of Sim to de- 
duce a new set of J^t values for four best known srf-shell 
emitters. The present estimation is based on the USD cal- 
culation with Vcoui plus Vq INC Hamiltonian, taken the 
average of the results obtained for various approaches to 
SRCs. The experimental ft values and the other correc- 
tions ((5roj S^ and (5ns) a-re taken from Ref. [5]. 

It is seen that for '^^Cl and ^''Ar, the deduced values 
are somewhat different from those adopted currently by 
Towner and Hardy. The implementation of our sd-shell 
results on Sjm is illustrated in Fig. 16, where corrected 
values of the 13 best known emitters from Ref. [5] are 
shown. Further analysis and calculation of the radial- 
overlap corrections is under way. 



B. Fermi Beta Decay to Non- Analogue States 

If the isospin-symmetry is broken, the Fermi beta decay 
between non-analogue 0^ states may take place. These 
transitions are of great interest since their rate is directly 
related to the degree of the isospin-mixing and it can help 
to test the model predictions [12]. Unfortunately, no ex- 
perimental data are known for the sd-shell nuclei. How- 
ever, we can explore the sensitivity of the matrix element 
of the non-analogue transitions (its isospin-mixing part) 
to the details of the shell-model Hamiltonian. In this 
context, we have performed calculations for two cases of 
^''Cl and '^^Ar, considered before by Ormand and Brown 
in Ref. [20] . The results are shown in Table XI, where we 
present calculations of Sjm for Fermi /^-transitions from 
the lowest 0+ state in "^^Cl or in ^"'Ar to the first excited 
0+ state in the respective daughter nuclei. Present calcu- 
lations are based on the Vcoui plus Vq INC Hamiltonian 
with strength parameters from Table V {Range 7). 

It is interesting to see that the values predicted for 
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■^^Cl by various interactions, including the early work of 
Ormand and Brown, are quite consistent, being roughly 
—2.7 3.1 X 10~^. However, the magnitudes of correc- 
tions predicted for a decay of '^'^Ar are much smaller and 
probably because of the sensitivity to the details of the 
calculations. Our result spread in a wide range between 
-8.8 X 10-5 and -48.2 x 10"^ for USD interaction and 
between -0.7 x 10"^ and -1.4 x lO'^ for USDA/USDB 
interactions. The value given in Ref. [20] is the small- 
est among the values, —3.5 x 10"^, while the parame- 
ters from Ref. [21] result contrary in the largest number, 
— 3.6x10-'^. It would be very interesting and useful to 
have experimental data to make a critical selection be- 
tween different predictions. 

VII. SUMMARY AND CONCLUSIONS 

In conclusion, we have presented a set of new empirical 
INC Hamiltonians in the sd shell-model space which accu- 
rately reproduces the isobaric mass splittings. The fitting 
procedure used in our work is close to that used earlier by 
Ormand and Brown, however, an advanced study of the 
harmonic oscillator parameters, modern approaches to 
SRCs, as well as, an updated and largerly extended exper- 
imental database and full sd-shell space calculations have 
been performed. In our model, besides a one-body term, 
an effective T = 1 component of the isospin-conserving 
interaction or the Wigncr term of the pion and p-meson 
exchange potentials have been exploited to model effec- 
tive charge-dependent forces of nuclear origin. More so- 
phisticated forms of those forces could be explored as well. 
The parameters of the INC part of the Hamiltonian were 
adjusted in a fit, designed to reproduce the known b and 
c coefficients of the IMME. Different types of the isospin- 
conserving interaction and various procedures to account 
for the SRC's lead to rather similar magnitudes of the 
rms deviations, with the best values around ~32 keV for 
rms of b coefficients and ^-9 kcV for rms of c coefficients. 
The quality of the fit and the Coulomb strength param- 
eters unambiguously suggest that in order to reproduce 
the experimental IMME coefficients, the electromagnetic 
interaction should be supplemented by nuclear charge- 
dependent forces. 

We believe that the constructed INC Hamiltonians can 
provide a high accuracy in the description of the isospin- 
symmetry forbidden processes. A few applications have 
been considered here with the purpose to demonstrate 
new features. First, we have been able to propose a quan- 
titative description of a staggering effect of the IMME 
b and c coefficients as a function of the mass number. 
This allowed us to conclude on the contribution of the 



Coulomb and nuclear charge-dependent forces to pairing. 

Second, we studied the validity of the IMME equation 
beyond the quadratic form in the lowest A = 32 quintet. 
Our calculations point towards the existence of a non-zero 
d coefficient in agreement with experimental data. The 
predicted values turn out to be very sensitive to the shell- 
model Hamiltonian used and thus more precise and more 
extensive data which may help to constrain theoretical 
parameters will be of great importance. 

Third, we present a new set of isospin-mixing correc- 
tions to sd-shell 0+ — 0+ beta decay rates. All Hamil- 
tonians provide surprisingly similar results, however, dif- 
ferent from the values of Towner and Hardy. A more 
advanced study of these corrections should be performed. 

Finally, preliminary calculations of the isospin-mixing 
corrections to the Fermi beta decay between non- 
analogue states suggest that these rates might be sensi- 
tive to the details of an INC Hamiltonian and thus could 
serve as a perfect test to a theoretical description. 

As a general conclusion, we hope that the Hamiltonian 
will be of large use and we intend to perform numerous 
applications to understand better predictions and differ- 
ences between various INC parametcrizations. Similar 
study has been performed for sdpf and pf model spaces, 
the results and various applications will be published else- 
where [94]. More extended and more precise experimen- 
tal data on isobaric masses, as well as measurements of 
the isospin-forbidden decay rates would be very helpful 
to test further and to refine our model. 
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TABLE XII. Comparison between fitted 6 coefficients witli experimental values of T = 1/2 doublets in sd-sliell nuclei. 



mass, A 


J"" 


6 (exp) 
(keV) 


b (fit)^ 
(kcV) 


References 


19 


1 + 
2 
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2 
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3 ^ 
2 
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3 + 
2 
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-t- 

2 
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■7 -t- 

2 
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23 


3 + 
2 
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5 ^ 
2 
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■7 -t- 

2 


4815.1 (11) 
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25 


5 + 
2 
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27 


2 
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1 

2 
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3 ^ 
2 
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1 

2 
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Q -t- 

2 


5835.0 (6) 


5760.3 




31 


1 + 
2 


6179.87 (98) 


6118.13 


[65, 99] 




3 
2 
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E 4- 

5 ^ 
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TABLE XIII. Comparison between fitted b- and c coefficients with experimental values of T = 1 triplets in sd-shell nuclei. 



mass, A 


J" 


(exp) 
(keV) 




b (fit)'' 
(keV) 


c (exp) 
(keV) 




c (fit)^ 
(keV) 


References 


18 


0+ 


3832.57 


(18) 


3842.20 


352.74 


(50) 


339.31 


[65, 95] 




2+ 


3785.2 


(2) 


3795.9 


267.1 


(5) 


268.4 




20 


2+ 


4216.2 


(5) 


4184.2 


185.1 


(20) 


193.7 


[65, 101] 




3+ 


4186 


(4) 


4184 


200 


(5) 


212 






4+ 


4206 


(4) 


4169 


181 


(5) 


186 




22 


0+ 


4594.74 


(16) 


4583.96 


312.03 


(26) 


297.64 


[65, 102] 




2+ 


4580.96 


(16) 


4577.16 


277.86 


(28) 


266.02 






4+ 


4570.0 


(2) 


4579.4 


230.6 


(4) 


230.0 






2+ 


4567.8 


(3) 


4581.0 


224.1 


(20) 


235.8 




24 


4+ 


4966.88 


(50) 


4956.56 


183.86 
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196.91 


[65, 97] 




1+ 


4943.7 
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4960.9 


182.0 
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2+ 
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0+ 
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3+ 
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[65, 99] 




2+ 


5670.57 


(55) 


5655.02 


179.61 


(56) 


186.36 




30 


0+ 


5967.64 


(150) 


5961.08 


275.93 


(153) 


272.91 


[65, 103] 




2+ 


5955.3 


(15) 


5944.0 


238.4 


(16) 


239.9 






2+ 


5919.7 


(15) 


5941.3 


220.7 


(16) 


229.8 




32 


1+ 


6267.2 


(5) 


6235.6 


194.1 


(6) 


199.8 


[65, 99] 




2+ 


6273.1 


(5) 


6232.9 


164.2 


(11) 


169.7 






0+ 


6243.9 


(5) 


6233.9 


149.2 


(11) 


153.2 




34 


0+ 


6559.48 


(17) 


6559.55 


285.50 


(18) 


271.69 


[65, 99] 




2+ 


6541.2 


(2) 


6542.3 


236.8 


(3) 


234.8 






2+ 


6551.1 


(3) 


6565.3 


198.1 


(4) 


204.3 






0+ 


6537 


(1) 


6574 


239 


(1) 


230 




36 


2+ 


6834.7 


(2) 


6859.3 


150.9 


(3) 


148.0 


[65, 99] 




3+ 


6845 


(1) 


6861 


224 


(1) 


220 






1+ 


6808.4 


(3) 


6854.1 


190.2 


(6) 


196.0 






1+ 


6843.5 


(3) 


6845.1 


239.9 


(7) 


235.8 




38 


0+ 


7110.5 


(2) 


7153.9 


283.7 


(4) 


268.6 


[65, 104] 




2+ 


7133.22 


(18) 


7102.99 


203.39 


(35) 


195.29 





^ {Vff^^)ijki of nuclear Hamiltonian USD had been used as the isospin-symmetry breaking term and UCOM SRC scheme was applied 
on Vcoul- 



[28] W. E. Ormand, Phys. Rev. C 53, 214 (1996). 
[29] W. E. Ormand, Phys. Rev. C 55, 2407 (1997). 
[30] B. J. Cole, Phys. Rev. C 54, 1240 (1996). 
[31] I. Hamamoto and H. Sagawa, Phys. Rev. C 48, R960 
(1993). 

[32] J. Dobaczewski and I. Hamamoto, Phys. Lett. B 345, 
181 (1995). 

[33] G. Colo, M. A. Nagarajan, P. V. Isacker, and A. Vitturi, 

Phys. Rev. C 52, Rll'75 (1995). 
[34] H. Sagawa, N. Van Giai, and T. Suzuki, Phys. Rev. C 

53, 2163 (1996). 
[35] H. Liang, N. Van Giai, and J. Meng, Phys. Rev. C 79, 

064316 (2009). 

[36] W. Satula, J. Dobaczewski, W. Nazarewicz, and 

M. Rafalski, Phys. Rev. Lett. 103, 012502 (2009). 
[37] W. Satula, J. Dobaczewski, W. Nazarewicz, and 



M. Rafalski, Phys. Rev. Lett. 106, 132502 (2011). 
[38] A. Petrovici, K. W. Schmid, O. Radu, and A. Faessler, 

Phys. Rev. C 78, 064311 (2008). 
[39] E. Caurier, P. Navratil, W. E. Ormand, and J. P. Vary, 

Phys. Rev. C 66, 024314 (2002). 
[40] N. Michel, W. Nazarewicz, and M. Ploszajczak, Phys. 

Rev. C 82, 044315 (2010). 
[41] N. Auerbach, Phys. Rev. C 79, 035502 (2009). 
[42] N. Auerbach, Phys. Rev. C 81, 067305 (2010). 
[43] E. Caurier, G. Martmez-Pinedo, F. Nowacki, A. Poves, 

and A. P. Zuker, Rev. Mod. Phys. 77, 427 (2005). 
[44] B. A. Brown and W. A. Richter, Phys. Rev. C 74, 

034315 (2006). 

[45] M. Honma, T. Otsuka, B. A. Brown, and T. Mizusaki, 

Phys. Rev. C 69, 034335 (2004). 
[46] F. Nowacki and A. Poves, Phys. Rev. C 79, 014310 



29 



TABLE XIV. Comparison between fitted 6 and c coefficients witli experimental values of T = 3/2 quartets in sd-shell nuclei. 



mass, A 


J'' 


6 (exp) 




b (fit)'' 


c (exp) 


c (fit)'' 


References 






(keV) 




(keV) 


(keV) 


(keV) 




19 


5 ^ 
2 


3982.7 


(40) 


3930.6 


240.7 (22) 


233.6 


[65, 95] 




Q -t- 

2 
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(50) 


3930.4 
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5 
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1 
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4399.6 


(23) 


4380.5 


216.8 (19) 
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ymmetry breaking term and UCOM SRC scheme was applied 


on Vcoul- 
















TABLE XV. 


Comparison 


between fitted b and 


c coefficients with 


experimental values of T 


— 2 quintets in 


sd-shell nuclei. 
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